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SECTION  1 
INTRODUCTION 

The  purpose  of  this  document  Is  to  provide  the  design  engineer  with  the 
Information  required  to  determine  the  heating  environments  and  certain  dynamic 
characteristics  of  sounding  rocket  payloads  r 5r  to  launch.  Depending  upon 
the  particular  mission,  there  are  certain  modifications  that  can  be  made  to 
an  existing  design.  These  modifications  could  possibly  eliminate  or  at  least 
considerably  reduce  certain  dynamics  and  heating  problems  associated  with  the 
reentry  payload.  The  ability  to  determine  these  potential  problems  prior  to 
launch  gives  the  engineer  a certain  amount  of  time  and  flexibility  In  consid- 
ering possible  defign  change  to  the  vehicle  to  prevent  these  problems  from 
occurring  during  flight. 

The  reentry  problem  that  must  be  addressed  by  the  design  engineer  Is 
complex.  Various  center  of  gravity  locations  combined  with  certain  release 
attitudes  can  often  result  In  complex  payload  attitude  dynamics  which  have  a 
pronounced  effect  on  the  reentry  heating  environment.  Because  of  these  vari- 
ations In  attitude,  constant  trim  methods  cannot  be  used  alone  to  adequately 
predict  the  trajectory  or  the  heating  environment  for  the  vehicle.  This  makes 
It  necessary  to  use  a combined  sophisticated  analysis  to  obtain  the  vehicle 
dynamics  and  heating  environment.  The  results  from  these  types  of  sophisti- 
cated methods  In  many  cases  are  too  complex  to  use  effectively  and  apply  only 
In  particular  situations.  A considerable  amount  of  effort  was  spent  to  assure 
that  the  results  presented  could  be  easily  used  and  general  enough  to  be  appli- 
cable to  numerous  reentry  heating  problems. 
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1.1  Scope  of  Contents 

The  results  presented  in  the  following  sections  are  applicable  to  the 
geometric  configurations  and  range  of  parameters  shown  in  Table  1.1.  The 
results  are  directly  applicable  only  to  cylinder  and  cone-cylinder  configura- 
tions with  a C.6.  offset  range  of  lOX  of  the  vehicle  length  forward  and  aft 
of  the  volumetric  C.G.  location.  The  initial  release  attitude  or  angle  of 
attack  range  is  from  90  to  180  degrees  referenced  to  the  earth's  local  ver- 
tical vector.  The  handbook  results  cover  a range  of  75  to  350  statute  miles 
in  release  altitude  and  consider  cylinder  diameters  from  22  to  38  Inches. 

The  skin  material  is  cons'  'ared  to  be  aluminum  with  skin  thicknesses  of  .125, 
.250  and  .500  inches.  Vehicle  fineness  ratios  of  2.5  to  5 are  considered 
for  the  Tjyllnder  configurations  while  fineness  ratios  of  3 to  7 are  consid- 
ered for  cone-cylinder  configurations.  The  resultant  ballistic  coefficient 
range  for  the  vehicles  is  between  10  and  250  1b/ft*.  The  results  presented 
In  the  following  sections  consider  neither  roll  nor  tumble  moments  Imparted  to 
the  vehicle  at  release. 


GEOMETRIC 

CYLINDER 

CONFIGURATION 

I 

C.  G.  Offset 
Fineness  Ratio 
Release  Attitude 
Release  Altitude 
Diameter 

Skin  Thickness 

Cone  Half  Angle 

-.10  < C.G*  < .1 0 

2.5  <£/d  < 5 

90“  <ai  < 180“ 

75  m1 . < Hi  < 360  mi . 
22  in.  < d < 38  in. 
.126  in.  < 5 < .500  in. 

-.10  < 

C.G.  ■ 

< .10 

3 < 

ih  < 

7 

90“  < 

(Hi  < 

180“ 

75  ml. 

< Hi 

< 350  ml. 

.22  in. 

< d < 

38  in. 

.125  1 

n.  < 6 

< .500  in. 

12“  < 

0 < 18 

0 

Table  1.1  Range  of  applicability  of  the  handbook  results 
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The  handbook  Is  comprised  of  essentially  eight  separate  sections  and  four 
appendices.  A concise  description  will  be  given  for  each  section  and  appendix: 

t Section  1,  In  addition  to  the  Scope  of  Contents  Section,  contains 
the  procedures  to  be  followed  when  solving  a problem  with  the  hand- 
book. The  solution  to  a sample  problem  Is  presented  and  can  be 
followed  as  an  example  when  solving  a reentry  for  an  actual  flight 
case. 

• Section  2 provides  all  the  necessa)7  information  required  to  de- 
termine the  trim  attitude  and  ballistic  coefficient  for  the  re- 
entry payload.  This  assumes,  of  course,  that  the  configuration 
type,  cross-sectional  area,  weight  and  C.G.  location  have  already 
been  specified  for  the  vehicle. 

• Section  3 contains  the  results  from  the  constant  trim  integrated 
heating  load  analysis.  This  heating  load  must  be  determined  in 
order  to  obtain  the  surface  temperatures  of  the  reentry  vehicle. 

• Section  4 presents  the  oscillation  and  tumble  correction  factors. 
These  factors  are  multiplied  by  the  constant  trim  heating  loads 
presented  1n  Section  3 to  correct  for  vehicle  oscillations  and 
tumbling  during  reentry, 

• Section  5 provides  all  the  necessary  Information  required  to  deter- 
mine the  surface  temperature  of  the  reentry  vehicle  once  the  tumble 
corrected  heating  load  lias  been  obtained. 

• Section  6 deals  specifically  with  the  dynamics  of  the  reentry  vehi- 
cle. Two  distinct  types  of  reentry  trajectories  are  considered, 
I.e.,  the  trajectories  wiilch  exhibit  normal  stable  flight  at  a 
constant  trim  attitude  and  those  which  consider  vehicle  oscillations 
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and  tumbling.  Methods  and  procedures  are  presented  1n  order  to 
determine  both  types  of  trajectories. 

f Section  7 presents  the  heating  rates  versus  altitude  and  velocity 
for  the  configuration  types  being  considered.  From  the  figures 
presented  In  Section  7,  the  heating  history  for  the  reentry  vehicle 
can  be  determined  once  the  trajectory  for  the  vehicle  has  been 
specified. 

• Section  8 presents  the  reentry  guidelines  that  should  be  followed 
when  using  the  handbook  to  solve  a vehicle  heating  or  dynamics 
problem, 

• Appendices  A-O  discuss  all  the  theoretical  methods  and  procedures 
used  In  determining  the  results  presented  In  the  following  sections. 

m 

1.2  Handbook  Usage  Procedures 

There  are  basically  three  types  of  analyses  which  can  be  performed  with 
the  handbook.  In  addition  to  the  heating  load  and  surface  temperature  analysis 
which  will  undoubtedly  be  the  most  Important  to  the  design  engineer,  trajec- 
tory and  heating  rate  analyses  can  also  be  performed.  The  basic  procedures 
Involved  In  using  the  handbook  arc  outlined  by  the  flowchart  presented  In 
Fig.  1.1 

A certain  amount  of  required  Information  Is  needed  to  begin  a reentry 
analysis.  This  Information  Includes: 

(1)  The  payload  geometric  configuration  type  and  fineness 
ratio 

(2)  Cone  nose  bluntness  ratio  If  cone-cylinder  geometry 

(3)  The  weight  of  the  payload 
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(4)  The  C.  G.  location  and  cylinder  cross-sectional  area 

4 

(5)  The  Initial  surface  temperature  of  the  vehicle  at  release 

^ (6)  The  Initial  attitude  and  altitude  of  the  payload  at  release 

(7)  Initial  roll,  tuntle  or  yaw  moments  Imparted  to  the  vehicle 
I at  release 

The  Initial  attitude  of  a reentry  vehicle  may  not  be  known  with  any  degree 
] . of  confidence  for  certain  missions.  In  these  situations,  engineering  Judgement 
and  knowledge  of  the  problem  must  be  used  to  come  up  with  the  best  estimate  of 

1 .1 

t;  this  attitude  since  It  Is  an  Important  factor  In  the  analysis. 

I i The  figures  In  Section  2 are  to  be  used  to  determine  the  trim  attitude  and 

S ! 

*V' 

'i  the  resultant  drag  coefficient  for  the  vehicle  type,  fineness  ratio  and  C.  G. 

I : 

I i location  given.  Once  the  drag  coefficient  has  been  established,  compute  the 

^ balHstll  coefficient  using  Eqn.  1.1. 

‘ Ballistic  Coefficient  “ W/CqA  1.1 

; where 

W “ weight  of  payload 

C[)  ■ drag  coefficient  based  on  cylinder  cross-sectional  area 

A “ cylinder  cross-sectional  area 

In  order  to  determine  the  surface  temperature  on  a vehicle,  the  integrated 
heating  load  must  be  obtained.  With  the  constant  trim  ballistic  coefficient, 
vehicle  trim  attitude,  and  release  altitude  already  established  at  this  point, 
determine  the  heating  load  from  Figs.  3.1  - 3.10  for  the  base  and  cylinder 
regions  of  the  vehicle.  Table  3.1  Is  provided  In  order  to  determine  the  heat- 
ing load  for  the  cone  section  when  considering  a cone-cylinder  geometry. 

Heating  load  magnitude  charts  arc  presented  in  Tables  3.2  and  3.3  for  the  con- 
figurations, C.  G.  locations,  and  release  attitudes  considered  In  the  6 DOF 
trajectory  and  tumble  heating  analyses.  The  heating  levels  presented  are  for 
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a release  altitude  of  200  miles.  While  the  magnitude  of  the  heating  loads 
presented  In  Tables  3.2  and  3.3  change  with  altitude,  the  distribution  does 
not.  From  Tables  3.2  and  3.3  the  engineer  can  determine  what  C.  G.  offset  and 
release  attitude  combinations  promote  the  highest  levels  of  heating  In  relation 
to  other  C.  G.  offsets  and  release  attitudes. 

The  heating  loads  obtained  In  Section  3 are  for  a constant  trim  reentry 
situation.  A tumble  correction  factor  must  be  applied  to  the  load  to  account 
for  vehicle  oscillations  and  tumbling.  Determine  this  tumble  correction  factor 
from  Figs.  4.1  -4.14  for  the  configuration  and  Initial  conditions  being  con- 
sidered. This  factor  is  multiplied  times  the  load  obtained  In  Section  3 to 
determine  a tumble  corrected  heating  load.  If  this  factor  Is  less  than  1,  It 
should  npt  be  applied  unless  the  Initial  release  attitude  is  known  with  a high 
degree  of  certainty. 

The  tumble  corrected  heating  loads  obtained  In  Section  4 will  be  used 
along  with  the  figures  presented  In  Section  5 to  determine  the  surface  tempera- 
ture of  the  reentry  vehicle.  Adjust  the  tumble  corrected  heating  load  for  sur- 
face temperature  variations  using  Fig.  5.1.  If  the  cylinder  diameter  for  the 
configuration  varies  from  24  Inches,  correct  the  heating  load  for  diameter 
change  using  F1g.  5.2.  If  the  vehicle  Is  rolling,  correct  the  stagnation 
line  heating  load  on  cylinder  section  by  the  roll  average  heat  transfer 
coefficient  using  F1g.  5.3.  Determine  the  vehicle  surface  temperature  using 
this  corrected  load  and  Fig.  5,4.  The  surface  temperature  obtained  from  Fig.  5.4 
considers  a constant  heat  capacity  for  the  vehicle  material.  Correct  this 
surface  temperature  for  variable  heat  capacity  using  Fig.  5.5. 

The  results  presented  In  Section  6 arc  used  to  determine  the  vehicle 
trajectory  once  the  ballistic  coefficient  has  been  established.  The  ballistic 
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coofficicdt  for  tho  mission  can  he  compared  l.o  t.hnse  presented  in  Table  G.l  in 
order  to  have  a check  on  the  ballistic  coefficient  calculation.  The  next  step 

in  determining  the  vehicle  trajectory  is  to  consider  the  results  in  Table  6.2. 

A tumbling  or  oscillating  trajectory  can  be  related  to  a constant  trim 
trajectory  using  these  results.  Determine  the  percent  of  variation  from  Table 
6.2  for  the  configuration  type  and  initial  conditions  being  considered. 
Calculate  an  equivalent  ballistic  coefficient  using  the  constant  trim  ballistic 
coefficient  and  this  percent  of  variation.  This  calculation  can  be  performed 
by  considering  the  following  relationship. 

/ V 

1.2 


W./w\.  V /w\ 
eq 


CT 


where 


/ W \ 
\l^f  eq 


The  equivalent  ballistic  coefficient  for  an  oscillating 
or  tumbling  trajectory 


f W 

I W CT 


Constant  trim  ballistic  coefficient 
V ■ Percent  of  variation  (from  Table  6.2) 

The  smaller  the  percent  of  variation  the  more  closely  the  trajectory  resembles 
a constant  trim  trajectory.  The  reentry  trajectory  of  the  vehicle  can  now  be 
determined  by  using  the  equivalent  ballistic  and  the  release  altitude  (from 
Figs.  6.1  - 6.9).  The  attitude  dynamics  for  the  reentry  vehicle  can  be 
established  for  the  particular  configuration  and  set  of  initial  conditions  by 
considering  the  information  provided  in  Table  6.3.  From  Table  6.3  it  is 
possible  to  determine  if  the  vehicle  is  oscillating  or  tumbling  or  whether  it 
is  trimmed  during  reentry.  It  is  also  possible  to  determine  if  these  oscilla- 
tions are  small  in  amplitude  or  whether  they  are  large  and  unstable.  Fig.  6.10 
has  been  provided  in  order  for  the  engineer  to  be  able  to  determine  the  maxi- 
mum descent  velocity  for  the  vehicle.  The  equivalent  ballistic  coefficient 
along  with  tha  release  attitude  can  be  used  in  Fig.  6.10  to  obtain  the  maximum 
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descent  velocity  of  the  reentry  vehicle  at  two  selected  recovery  altitudes. 

Section  7 deals  specifically  with  the  heating  rate  environments  for  the 
reentry  payloads.  If  a heating  rate  analysis  1s  being  considered,  examine  the 
magnitude  of  the  percent  of  variation  determined  In  Section  6.  If  the  absolute 
value  of  this  percentage  Is  greater  than  5X,  this  Is  an  Indication  that  vehicle 
oscillations  or  tumbling  will  occur.  The  heating  rates  presented  In  Figs.  7.1  - 
7.10  are  for  constant  trim  attitudes.  The  heating  rate  history  for  a vehicle 
cannot  be  determined  using  these  figures  If  the  vehicle  attitude  Is  not  known 
at  each  specific  time  point  which  Is  the  case  1n  an  oscillating  or  tumbling 
trajectory.  If  the  percent  of  variation  Is  larger  than  5X  and  the  heating  rate 
history  of  the  vehicle  Is  required,  a complete  6 OOF  trajectory  analysis  and  a 
tumbling  heating  analysis  must  be  performed  for  these  results.  If,  however, 
the  absolute  value  of  the  percentage  Is  less  than  591,  this  means  that  the 
vehicle  Is  essentially  trimmed  at  a constant  attitude.  For  these  types  of 
situations,  select  a heating  rate  plot  from  Figs.  7.1  - 7.10  which  Is  closest 
to  the  trim  attitude  of  the  vehicle  and  superimpose  the  reentry  trajectory 
that  was  determined  In  Section  6.  The  heating  rate  history  can  then  be 
determined  for  both  the  cylinder  and  base  regions  of  the  payload.  Table  7.1  Is 
provided  In  order  to  determine  heating  rates  on  the  cone  section  for  the  cone- 
cylinder  configurations.  Fig.  7.11  is  also  presented  In  order  to  show  the 
altitude  cutoff  limits  for  the  heating  rate  analysis.  These  are  the  altitudes 
at  which  aerodynamic  heating  ends  and  cooling  begins. 

The  reentry  guidelines  are  presented  In  Section  8.  When  potential  heating 
or  dynamic  problems  are  Indicated  for  a particular  flight  case,  follow  the 
recommendations  presented  in  Section  8 prior  to  changing  the  vehicle  design. 
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1.3  Sample  Problem 

A detailed  step  by  step  solution  to  a sample  problem  Is  presented  In  this 
section.  The  problem  Illustrates  the  techniques  and  procedures  that  a user 
should  follow  In  solving  a problai.  This  problem  Is  representative  of  the 
types  of  problems  which  will  most  likely  be  encountered  while  performing  a 
flight  reentry  heating  analysis. 

PROBLEM: 

Determine  the  maximum  descent  velocity  at  60K  feet,  the  attitude  dynamics, 
maximum  heating  rate  and  surface  temperature  on  the  cylinder  region  of  a blunt 
cone-cylinder  configuration. 


CONFIGURATION  AND  MISSION  CHARACTERISTICS: 


Item 

Symbol 

Value 

Units 

Configuration  type 

cone-cylinder 

Cylinder  diameter 

d 

37.5 

in. 

Vehicle  length 

« 

158 

in. 

Vehicle  weight 

w 

1742 

lbs . 

Nose  bluntness  ratio 

0.132 

Skin  thickness 

0.250 

in. 

Release  altitude 

1 75 

statute  miles 

C.  G.  location 

x/l 

.645 

Release  attitude 

«1 

90 

degrees 

cone  half  angle 

0 

14 

degrees 

Vehicle  roll  rate 

r 

- k “*1 

.5 

rps 
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SOLUTION  FOR  SURFACE  TEMPERATURE: 


Step  No. 


1 

Calculate  fineness  ratio 

i!/d  “ 4.21 

2 

From  Fig.  2.2  determine  volumetric  C.  G. 
location 

x/i  = .620 

3 

Calculate  vehicle  C.  G.  offset  from 
volumetric  C.  G. 

C.  G.  Offset  = .025  l 
AFT  of  Vol.  C.G. 

4 

From  Fig.  2.3  determine  planform  area 
ratio 

Handbook  results  are 
valid  (Section  2) 

5 

From  Fig.  2.4  determine  trim  attitude 

ay  - lei" 

6 

From  Fig.  2.5  determine  drag  coefficient 

Cl)  - 1.9 

7 

Calculate  cylinder  cross-sectional  area 

A - 7.67  ft^ 

8 

_ Calculate  ballistic  coefficient 

W/c^A  " 115.5  Ibm/ft^ 

9 

From  Fig.  3.8  and  3.9  linearly  inter- 
polate for  heating  load 

Q^.g5  « 207  Btu/ft2 

10 

Calculate  attitude  offset  parameter 

P " -.441 

*11 

Calculate  vehicle  fineness  ratio  as  if 
the  vehicle  were  a sharp  cone-cylinder 
geometry,  ia  “ 0 

i/d  “ 4.48 

12 

Determine  from  Fig.  4.6  and  4.7  tumble 
correction  factor  (linearly  Interpolate 
with  vehicle  fineness  ratio  from 

Step  11) 

Qt/Qct  - 1*87 

13 

Calculate  tumble  corrected  heating  load 

387  Btu/ft2 

14 

Determine  from  Fig.  5.1  variable  surface 
temperature  corrected  load 

Q “ 360  Dtu/ft2 

15 

Determine  diameter  adjustment  factor 
from  Fig.  5.2 

Factor  = .91b 

16 

Adjust  load  from  Step  14  for  diameter 
change 

Q » 329  Btu/ft2 

17 

Determine  the  factor  Fr/h.n  from  Fig. 

5.3  for  ax  and  multiply  times  Qi»o 

Peso  = 145 

18 

Determine  constant  heat  capacity  surface 
temperature  from  Fig.  5.4 

Tsc  “ 843»R 

n 
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solution  for  surface  temperature  (Continued): 
Step  No. 


#*19 

Correct  surface  temperature  for 
variable  heat  capacity  using  Fig.  5.S 

Ts  “ 833*8 

SOLUTION  FOR  MAXIMUM  VELOCITY  AT  60K  FEET: 

Step  No. 

20 

Determine  percent  of  variation  from 

Table  6.2  (linearly  Interpolate  between 
fineness  ratios  using  fineness  ratio 
from  Step  1 1 ) 

Percent  » -7.8* 

21 

Calculate  equivalent  ballistic  coeffi- 
cient (use  Eq.  1.2) 

(W/cqA)«‘1  “110 

22 

From  Fig.  6.4  determine  maximum 
velocity  at  60K  feet 

Maximum  velocity 
« 4000  ft/sec 

• 

SOLUTION  FOR  ATTITUDE  DYNAMICS: 

Step  No. 

23 

From  Table  6,3  determine  attitude 
dynamics 

Vehicle  oscillations 
are  less  than  25®  about 
Of.  Vehicle  Is  re- 
entering aft  end  first 

SOLUTION 

FOR  MAXIMUM  HEATING  RATE: 

Step  No. 

24 

The  percent  of  variation  from  Step  19  Is 
7.8*  which  Is  greater  than  the  5*  limit 
Imposed.  An  approximate  maximum  heating 
rate  will  try  to  be  obtained  even 
though  the  vehicle  Is  not  at  constant 
trim 

25 

Superimpose  the  trajectories  from  Fig. 

6.4  onto  the  heating  rate  plots  for  the 
cylinder  section  In  F1g.  7.8 

26 

Determine  maximum  heating  rate 

Maximum  heating  rate 
“ 20  Btu/ft^-sec 

12 
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♦This  calculation  must  be  performed  since  the  results  in  Section  4 were 
determined  for  sharp  cone-cylinder  geometries  and  blunt  cylinder  geometries 
only.  Determine  the  length  of  the  cone-section  from  Its  theoretical  apex 
(see  Fig.  2.1  and  Eq.  2.1). 

‘c  ■ 

Calculate  cone  length 

la  - .132  ig*  g.926  In. 

Determine  vehicle  length  from  theoretical  apex 

1 ♦ ij  ■ 168  + 9.926  « 167.926  In. 

Calculate  fineness  ratio  from  theoretical  apex 

f - 167.926/37.5  » 4.48 

♦♦This  Is  the  final  estimate  of  the  surface  temperature  on  the  cylinder  section. 
This  value  assumes  an  Initial  surface  temperature  at  apogee  of  660''R. 
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SECTION  2 

AERODYNAMIC  CHARACTERISTICS 

The  data  based  theoretical  methods  discussed  in  Appendix  A provide  the 
necessary  aerodynamic  coefficients  used  In  both  the  six-degree-of-freedotn  tra- 
jectory simulation  analysis  and  the  three-degree-of-freedom  ballistic  reentry 
analysis.  The  normal  force,  axial  force,  and  pitching  moment  coefficients  were 
computed  using  these  methodologies.  The  damping  coefficients  which  were  an 
Integral  part  of  the  6 DOF  trajectory  analysis  were  computed  using  Ref.  1.  The 
aerodynamics  computed  using  these  theoretical  methods  consider  a sharp  cone 
geometry  for  the  cone  section  of  the  cone-cylinder  configurations,  The  sharp 
cone  georffetry  considered  In  the  handbook  analysis  can  be  defined  as  one  In 
which  the  cone  tip  and  theoretical  apex  coincide.  While  It  Is  true  that  the 
aerodynamics  are  not  directly  applicable  to  blunt  cone-cylinder  geometries, 
there  are  certain  degrees  of  nose  bluntness  for  which  the  handbook  results  are 
valid.  In  order  to  define  nose  bluntness  consider  the  schematic  In  Fig.  2.1. 


Fig.  2.1  Cone-cylinder  configuration  with  cone  nose  bluntness 
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Tho  nosi;  bluntness  rnUo  will  Lio  ilorinud  ds  t.lu'  d i stonco  rroiii  tho  thoorotical 
apox  to  tho  di;t.iidl  noso  r.itiooii  to  tho  distanco  from  tho  thoorotlcfll  apox 
to  tho  i;on(!  Iiaso. 

. ^ . ^atanl  2.1 

where 

“ distance  from  the  theoretical  apex  of  the  cone  to  the  cone  base 

■ distance  from  the  theoretical  apex  to  tho  hemisphere-flat  face 
tangent  point 

In  order  to  maintain  consistency  In  tho  application  of  the  handbook  results  to 
both  sharp  and  blunt  cone  cylinder  geoinotrios,  tho  vehicle  fineness  ratio  will 
always  be  defined  unless  otherwise  stated  as  tho  distance  from  the  nose  of  the 
vehicle  to  tho  base  of  the  cylinder  afterbody. 

A volumetric  C.  G.  location  must  be  determined  for  tho  reentry  vehicle  as 
a first  step  In  the  reentry  analysis  regardless  of  the  vehicle  geometry.  Fig. 
2.2  presents  the  volumetric  C.  G.  locations  for  cone-cylinder  geometries  with 
nose  Ijluntnoss  ratios  ranging  from  zero  for  a sharp  conc-cylinder  geometry  to 
one  for  a cylinder  geometry.  When  a blunt  nose  cone  geometry  Is  being  con- 
sidered In  a reentry  analysis,  tho  nose  bluntness  ratio  Is  first  determined  by 
using  oguatlon  ZA,  The  volumetric  C.  0.  location  can  then  be  deternilnod  for 
tho  vehicle  fineness  ratio  and  specific  nose  bluntnoss  ratio  using  Fig.  2.2. 

As  the  COMO  ttose  bluntnoss  increases  the  vetilcle  aerodynamics  begin  to 
cliange.  There  are  certain  nose  bluntness  cutoff  limits  that  must  bo  considered, 
lioyohd  those  limits  tlie  velilclo  aerodynamics  liave  changed  to  the  extent  that 
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the  handbook  results  are  no  longer  valid  for  the  particular  geometry.  When 
considering  a blunt  cone-cylinder  configuration,  the  bluntness  ratio  determines 
whether  the  aerodynamics  are  more  closely  associated  with  those  of  a sharp  cone 
geometry  or  those  of  a cylinder  geometry.  Therefore,  basically  tv/o  cutoff 
limits  exist.  One  limit  would  be  the  bluntness  ratio  that  would  render  the 
sharp  cone-cylinder  results  invalid,  and  the  other  limit  would  be  that  blunt- 
ness which  would  make  the  cylinder  results  invalid.  These  limits  will  change, 
of  course,  with  vehicle  fineness  ratio.  It  is  not  possible  to  determine  these 
exact  limits  without  performing  an  in-depth  analysis.  It  is  possible,  however, 
to  make  some  generalized  approximate  statements  regarding  these  limits.  The 
aerodynamics  for  a vehicle  are  direct  functions  of  the  planform  area  of  that 
vehicle.  ^As  the  planform  area  changes  there  are  corresponding  changes  in  the 
aerodynamics.  Fig.  2.3  presents  the  planform  area  for  a blunt  cone-cylinder 


geometry  ratioed  to  the  planform  area  of  a sharp  cone-cylinder  geometry.  This 


ratio  (Ahr/Acr)  is  presented  versus  fineness  ratio  and  nose  cone  bluntness 


ratio.  The  fineness  ratio  for  this  specific  case  will  be  the  distance 


from  the  theoretical  cone  apex  to  the  base  of  the  cylinder  afterbody.  When  con- 


sidering a blunt  nose  cone-cylinder  geometry  for  a specific  flight  case,  it  is 
generally  felt  that  the  handbook  results  should  be  valid  and  accurate  if  the 
planform  area  ratio  is  greater  than  0.95.  The  data  book  will  become  less  and 
less  accurate  for  planform  area  ratios  less  than  0.95.  A generalized  statement 
regarding  a practical  cutoff  limit  for  a cylinder  configuration  cannot  be  made 
as  easily  however.  It  has  been  experimentally  determined  and  theoretically 


verified  that  the  addition  of  a hemispherical  cap  to  one  end  of  a cylinder  can 


greatly  alter  the  aerodynamics  of  the  cylinder.  For  this  reason,  it  is  recom- 
mended that  blunt  cylinder  configurations  be  separately  analyzed  as  to  their 
reentry  characteristics. 
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A trim  attitude  Is  established  for  a reentry  vehicle  when  Its  C.  G.  location 
and  aerodynamic  force  center  conclde.  This  1s  the  point  at  which  there  are  zero 
moments  and  torques  on  the  reentry  body.  The  aerodynamic  force  center  versus 
angle  of  attack  can  be  computed  using  the  following  equation. 


where 

a Aerodynamic  force  center  location  non-dimensional Ized  by  the 
vehicle  length 

Vm 

^ « Pitching  moment  reference  point  location  non-dimensional Ized  by 

the  vehicle  diameter 

C|„  = Pitching  moment  coefficient 

C[^  • Normal  force  coefficient 

^ ® Vehicle  diameter  to  length  ratio 

Aerodynamic  force  centers  were  computed  using  Equation  2.2  for  the  configuration 
types  and  range  of  fineness  ratios  presented  In  Table  1.1.  An  angle  of  attack 
range  of  5 to  175  degrees  was  considered  In  the  calculations.  The  aerodynamic 
force  center  Is  undefined  at  0°  and  180°  when  Equation  2,2  Is  used.  This  Is 
d le  to  the  fact  that  both  the  normal  force  and  pitching  moment  coefficients  are 
zero  at  these  angles  of  attack.  The  aerodynamic  force  centers  versus  fineness 
ratio  for  constant  trim  attitudes  are  presented  In  Fig.  2.4  for  the  cone- 
cylinder  configurations.  The  same  results  are  presented  In  Fig.  2.6  for  the 
cylinder  configurations.  Aerodynamic  drag  coefficients  for  the  cone-cylinder 
configuration  versus  trim  attitude  and  fineness  ratio  are  given  In  Fig.  2.5 
while  similar  results  arc  presented  for  the  cylinder  configuration  Is  Fig,  2.7. 
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C,G.  Location,  X/t 


Trim  Attitude,  Oj  (degrees) 

coefficient  versus  trim  attitude  and  fineness  ratio  for  cylinder  configurations. 


SECTION  3 

INTEGRATED  HEATING  LOADS 

The  heating  loads  discussed  in  this  section  fall  basically  into  two  cate- 
gories. Ballistici  or  constant  trim  entry  heating  loads,  are  discussed  in  Sec- 
tion 3.1.  Relative  heating  load  magnitudes  for  a number  of  trajectories  which 
include,  not  only  constant  trim  entries,  but  also  oscillating  and  tunbling 
type  entries  are  discussed  in  Section  3.2. 

3,1  Constant  Trim  Entry  Loads 

Integrated  heating  loads  were  computed  using  the  Bentry  Computer  Program 
(Ref.  C.1  - Appendix  C)  for  each  of  the  constant  trim  ballistic  trajectories 
presented  in  Figs.  6.1  - 6.9  of  Section  6.  The  heating  loads  were  computed  for 
a release  altitude  range  of  50  to  350  statute  miles.  The  results,  presented 
in  Figs.  3.1  - 3.10,  are  applicable  to  the  base  and  cylindrical  sections  of 
both  the  cylinder  and  cone-cylinder  configurations. 

The  heating  load  results  presented  in  Figs.  3.1  to  3.10  were  calculated 

for  a reference  diameter  of  24  inches,  for  angles  of  attack  of  90  to  180  de- 

grees using  a reference  temperature  of  660“R.*  The  heating  load  may  be  adjusted 
for  diameter  variations  using  Fig.  5.2.  To  compute  a load  for  trim  angles  be- 
tween 0 and  90  degrees  use  the  following  relations: 

Cylinder:  a-j-  ■ 180  - a 

Base:  * 90 

A wall  temperature  adjustment  to  the  heating  load  is  provided  in  Section  7 by 
Fig.  7.1. 

*660*’R  was  assumed  to  be  a nominal  estimate  of  the  wall  temperature  at  apogee. 
Differences  in  wall  temperatures  at  apogee  are  discussed  in  Section  8. 
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The  heating  load  to  the  cone  section  of  a cone-cylinder  body  may  be  calcu- 
lated  using  the  algorithm  given  In  Table  3.1  where  the  cylinder  Is  used  as  a 
reference. 

I 

3.2  Oscillating  and  Tumbling  Entry  Heating  Levels 

'4  ' 

The  Bentry  program  has  the  capability  of  determining  heating  rates  and 
integrated  loads  on  tumbling  reentry  bodies  once  the  trajectory  for  the  tum- 
bling body  has  been  established.  The  matrix  of  runs  considered  in  the  6 DOF 
trajectory  analysis  Is  presented  in  Table  6.2  of  Section  6.  Integrated  heat- 
ing loads  were  computed  using  the  Bentry  program  for  each  of  these  trajectories 
which  Included  those  with  vehicle  oscillations  and  tumbl1ng»  as  well  as  the 
constant  trim  reentries.  While  these  Bentry  runs  for  heating  loads  were  per- 
formed nfSInly  for  determining  the  oscillating  and  tumbling  correction  factors 
presented  in  Section  4,  they  are  also  used  In  determining  relative  heating  load 
magnitudes  for  various  combinations  of  C.G.  offsets  and  Initial  release  atti- 
tudes for  the  payload  configurations  considered  In  the  analysis.  These  results 
are  shown  In  the  heating  load  magnitude  charts  presented  In  Tables  3.2  and  3.3 
A number  of  trajectories  were  determined  with  the  6 DOF  analysis  for  Initial 
release  altitudes  other  than  200  miles.  From  these  trajectory  runs,  It  was 
determined  that  the  vehicle  dynamics  with  respect  to  the  attitude  of  the  ve- 
hicle during  an  oscillating  or  tumbling  reentry  were  essentially  constant  for 
the  release  altitudes  examined.  The  extensive  trajectory  analysis  1n  which  C.G. 
locations  and  release  attitudes  were  varied  and  for  which  the  heating  load  cal- 
culations were  performed  considered  a 200  mile  release  altitude.  Therefore, 
only  the  overall  magnitude  of  the  heating  loads  presented  In  Tables  3.2  and  3.3 
would  change  with  release  altitude  and  not  the  distribution  since  the  vehicle 
attitude  dynamics  are  essentially  constant  with  release  altitude. 
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A number  of  trends  can  be  established  and  a number  of  conclusions  made 
regarding  the  results  presented  In  Tables  3.2  and  3.3.  it  is  obvious  from 
Table  3.2  that  the  base  load  environment  Is  considerably  higher  In  practically 
all  cases  where  the  C.G.  offset  Is  aft  of  the  volumetric  C.G.  location.  This 
occurs  since  the  base  region  is  subjected  to  the  flow  In  situations  where  the 
vehicle  Is  trimmed  or  Is  oscillating  aft  end  first  during  reentry,  Notice 
that  the  maximum  loads  on  the  base  occur  for  the  4,  5 and  7 fineness  ratio 
cone-cylinder  configurations  at  C.G.  offsets  of  2.55!  forward,  2.55!  aft  and  55! 
aft,  respectively.  These  maximum  levels  are  essentially  Independent  of  Initial 
release  attitude.  For  the  cylinder  section  of  the  cone-cylinder  configuration, 
a totally  different  heating  load  level  distribution  is  observed.  The  1055  for- 
ward location  gives  the  highest  load  levels  and  the  distribution  Is  highly 
dependent  on  the  Initial  release  attitude,  but  essentially  Independent  of  fine- 
ness ratio.  For  the  base  and  cylinder  regions  of  the  cylinder  configurations, 
the  maximum  load  levels  occur  for  a 105!  C.G.  aft  offset.  Notice  the  difference 
In  the  load  levels  on  the  base  region  between  the  cone-cylinder  configurations 
of  fineness  ratios  4 and  5 and  the  same  fineness  ratio  cylinder  configurations. 
Maximum  heating  loads  occur  on  the  base  region  of  both  types  of  configurations 
when  the  vehicle  Is  trimmed  aft  at  140®  angle  of  attack.  This  Is  not  necessar- 
ily the  case,  however,  when  a vehicle  Is  oscillating  about  a 140®  trim  point 
since  the  base  may  be  subjected  to  flow  more  often  at  angles  of  attack  which 
vary  from  140®.  These  variations  are  dependent  upon  the  amplitude  of  the  oscil- 
lations. It  Is  entirely  possible  during  an  oscillating  trajectory  that  oscil- 
lations about  a 130®  or  150®  trim  attitude,  for  Instance,  could  give  the  highest 
heating  levels  since  the  base  would  then  be  subjected  more  often  during  the  tra- 
jectory to  the  140®  maximum  load  attitude.  For  the  cone-cylinder  configurations 

with  fineness  ratios  of  4 and  5,  with  C.G.  offsets  of  2.55!  forward  and  2.5;t  aft, 
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respectively,  the  constant  trim  attitudes  (from  Fiqs.  2.2  and  2.4  in  Section  2) 
are  147*  and  158®.  These  C.G.  offsets  result  in  the  highest  load  levels  for 
the  4 and  5 fineness  ratio  cone-cylinder  configurations.  Similarly,  constant 
trim  attitudes  of  132®  and  138®  for  10%  aft  offsets  in  C.R.  give  the  largest 
loads  for  the  4 and  5 fineness  ratio  cylinder  configurations. 

A number  of  other  comparisons  can  be  made  and  trends  established  concern- 
ing relative  heating  levels  with  Tables  3.2  and  3.3  for  each  of  the  cases  and 
configurations  considered  in  the  analysis.  The  results  presented  in  Tables 
3.2  and  3.3  should  prove  to  be  valuable  to  the  engineer  in  that  he  can  deter- 
mine very  quickly  and  without  any  calculations  the  relative  magnitude  of  the 
heating  levels  anticipated  for  his  particular  mission. 
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TABLE  3.1 

CONE  HEATING  LOADS 


The  procedure  for  computing  heating  loads  to  the  cone  of  a cone-cylinder 
configuration  is  as  follows: 

(1)  Determine  the  angle  of  attack  of  the  cone's  windward  surface 

a. 


(2) 

(3) 


cone 


cyl  c 


a , ■ Angle  of  attack  of  cylinder 

cyl 

» Semi -vertex  angle  of  the  cone 


^^«cone>  ^80“.a,,„3 


180' 


Determine  the  heating  load  for  a cylinder  at  iv  from  Figs.  3.1 

cone 

to  3.10. 

Determine  the  cone  radius  at  the  mid  length  of  the  cone. 


(4)  Determine  the  load  adjustment  factor  for  the  equivalent  cone  dia- 
meter, 2r  , by  using  Fig.  B.?. 

C 

(5)  Make  additional  adjustments  as  required  for  a cylinder  calculation. 


Cylinder  Heating  Load  ^bg^1^e  Chart  (H  = 200  rai. 


1000 


tim 


Fig.  3,2  Heating  Loads  for  a^=100' 
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Fig.  3-8  Heating  Loads  for  a-=]60‘ 


Heating  Loads  for  a_=180' 
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SI.CT10N  4 

OSCILLATlNd  AND  lUMDI.I.  CORRI.CTIDN  I'ACTORS 

As  discussed  in  previous  sections,  various  C.  (5.  locations  and  initial 
release  attitude  combinations  cause  large  vehicle  oscillations  and  tuniblino 
during  reentry  for  tl)e  types  of  geometric  configinwtions  being  considered  in 
this  analysis.  The  effects  of  those  types  of  dynamics  on  tlie  vehicle's  entry 
tf'ajectory  as  to  its  altitude,  velocity,  and  time  (listory  were  discussed  in 
Section  6.  The  aerodynamic  heating  on  a reentry  vehicle  is  determined  not 
only  by  the  vehicle's  velocity  and  altitude,  but  also  its  attitude.  Although 
the  variations  in  the  vehicle's  descent  velocity  cliungos  the  licating  load,  the 
variation^  in  attitude  have  a more  pronounced  effect.  Ttiis  can  be  seen  by 
examifilmj  the  heating  load  plots  given  In  Figs.  3.1  - 3.10  in  Section  5 for 
the  range  of  90  to  IHU  degrees  in  ri'ontry  angle  of  attack. 

During  typical  large  oscillating  trajeclorii's , tin;  anglo  of  attack  may 
vary  180"  for  one  complete  oscillation  of  the  reentry  vehicle,  This  variation 
could  be  3()0‘’  when  tlu'  vehicle  completely  tumbles.  The  lieating  rate  on  thi' 
reentry  vehicle  can  vary  from  a very  small  rate  to  a tremendously  large  rate; 
during  one  complete  oscillation.  Therefore,  in  ordi'r  to  adecpiately  pei'foriu  a 
heating  analysis,  the  inteejrated  heating  load  must  be  determined  which  consider 
this  large  variation  in  vehicle  attitude.  This  is  the  rationale  in  usin(| 
tumbling  entry  heating  program  such  as  the  IRnilry  Compuier  Program  for  the 
heating  analysis. 

Tables  3.1  and  3,2  in  Section  3 present  the  n^lative  levels  of  lUMting  for 
the  cent j (luration  tyiuis  and  C.  0.  offset  and  attitude  combinations  considered 
in  tills  analysis,  file  tra jeetories . for  whicli  thes('  loads  weri'  nlitiiined,  cover 
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practically  all  possible  types  of  constant  trim,  oscillatinq,  and  tumbling 
reentry  trajectories.  The  heating  loads  computed  with  the  Bentry  Program  for 
these  trajectories  were  non-dimensional ized  by  the  loads  obtained  for  a constant 
trim  ballistic  reentry.  Thesu  non-dinujnslonal ized  values  which  are  called 
tumble  correction  factors  were  plotted  versus  attitude  offset  parameter,  P,  and 
presented  in  Figs.  4.1  - 4.14.  The  attitude  offset  parameter  is  defined  by 
equation  4.1 

' (It  . a 

p.-i-I  4.1 

a-j- 

where 

a.j  = The  initial  vehicle  release  attitude  (Deg.) 

uj  “ The  vehicle  trim  attitude  determined  from  Section  2 (Deg.) 

It  is  obvious  from  equation  4.1  that  P Is  equal  to  0 if  the  vehicle  were  released 
at  its  trim  attitude  and  that  the  value  of  P becomes  larger  the  further  the 
release  attitude  is  from  the  vehicle  trim  attitude. 

It  is  Important,  in  order  to  correctly  apply  the  tumble  correction  factors, 
that  t.he  philosophy  behind  the  development  of  these  factors  bo  discussed.  A 
trim  attitude  can  be  established  for  a reentry  vehicle  by  the  methods  presented 
in  Section  2 if  the  aerodynamic  force  center  and  C.  G.  location  is  known  for  the 
vehicle.  If  the  vehicle  trims,  it  will  always  trim  at  this  attitude  and  If  the 
veliicle  oscillates  or  tumbles,  it  will  always  oscillate  or  tumble  about  this 
trim  attitude.  Therefore,  regardless  of  what  type  of  vehicle  attitude  dynamics 
are  involved,  a reference  heating  load  can  be  computed  as  if  the  vehicle  were 
triinmed  at  its  theoretical  trim  attitude.  This  is  precisely  the  manner  in  which 
the  constant  trim  loads  were  computed  in  order  to  non-dimensional i ze  the  heating 
loads  computed  with  the  Bentry  Program  and  to  form  the  tumble  correction  factors 
presented  in  this  section. 
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The  tumble  correction  factors  presented  in  Flqs.  4.1  - 4.14  were  deter- 
mined for  all  the  cases  and  configurations  considered  In  this  analysis.  If 
the  fineness  ratio,  C.  G.  location,  or  initial  release  attitude  for  a flight 
case  varies  from  the  cases  considered,  linear  interpolate  in  order  to  determine 
the  tumble  correction  factors  for  that  case. 

Tumble  correction  factors  of  20  are  not  uncommon  for  a number  of  cases 
considered.  Even  though  these  factors  are  large,  the  resultant  tumble  corrected 
heating  loads  obtained  when  using  these  factors  are  small.  This  is  due  to  the 
fact  that  the  constant  trim  load,  is  very  small.  For  instance,  notice 
in  Fig.  4.1  that  a factor  of  21  is  obtained  for  the  base  region  of  a 7 fineness 
ratio  cone-cylinder  configuration  with  a 2,B%  offset  forward  in  C.  G.  location 
released  at  an  attitude  of  180®.  The  constant  trim  load  for  this  case  is 
9 Btu/ft^.  The  tumble  corrected  load  is  then  equal  to  this  load  times  the  tumble 
correction  factor  of  21  which  is  equal  to  189  Btu/ft^  which  is  a comparatively 
small  heating  load.  This  is  usually  the  case  when  large  factors  are  obtained. 

Even  though  the  tumble  correction  factors  are  presented  for  use  in  the  section, 
one  should  exercise  great  care  in  their  application.  This  is  especially  true  in 
applying  tumble  correction  factors  less  than  one.  For  engineering  design  purposes, 
the  following  gu 1 delines  are  recommended.  Apply  the  correction  factors  only  when 
they  are  greater  than  one.  Use  a factor  of  one  in  situations  where  the  tumble 
correction  factor  is  less  than  one.  In  this  way  a certain  amount  of  conservatism 
is  maintained  in  the  results  since  the  heating  load  is  never  reduced  to  a value 
less  than  the  constant  trim  heating  load. 
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SECTION  5 

SURFACE  TEMPERATURE  CALCULATIONS 

A considerable  amount  of  discussion  concerning  reentry  trajectories 
and  heating  rates,  and  loads  on  reentry  vehicles  has  been  presented  In  the 

preceding  sections.  This  section  deals  with  the  methods  and  procedures 

1 

Involved  1n  using  the  heating  loads  determined  In  these  sections  to  obtain 
vehicle  surface  temperatures  during  reentry.  The  maximum  vehicle  surface 
temperatures  are  of  main  concern  to  the  design  engineer.  From  these  tempera- 
tures he  can  determine  whether  the  vehicle  structure  and  skin  can  withstand 
the  reentry  environment.  He  could  also  determine  If  the  surface  temperatures 
would  be<of  great  enough  magnitude  that  the  vehicle  could  not  be  safely  reused 
In  future  missions. 

As  discussed  In  Section  3,  and  more  thoroughly  discussed  1n  Appendix  C, 
the  Integrated  loads  presented  In  Figs.  3.1  - 3.10  consider  a reference 
surface  temperature  of  660*R  In  the  load  calculation  procedures.  In  reality, 
the  actual  heating  load  varies  with  vehicle  surface  temperature.  Since  the 
surface  temperature  Is  a function  of  the  vehicle  skin  thickness,  the  heating 
load  becomes  a function  of  skin  thickness.  Fig.  5.1  presents  the  actual 
heating  load  (Q)  versus  the  heating  load  computed  at  the  reference  temperature 
of  660'’R  (Qeeo)  fo**  wall  thicknesses  of  .125,  .250  and  .500  Inches.  After 
a heating  load  1$  obtained  In  Section  3 and  then  corrected  for  a possible 
oscillating  or  tumbling  reentry  In  Section  4,  the  variable  surface  temperature 
correction  can  be  made  and  the  actual  load  determined  by  using  Fig,  5.1. 
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The  aerodynamic  heating  rates  and  subsequent  loads  on  reentry  vehicles 


are  dependent  on  the  vehicle  radius  or  diameter.  In  general,  the  heating 
analysis  was  performed  for  a cylinder  diameter  of  24  Inches  for  both  the 


cylinder  and  cone-cylinder  configurations.  Both  the  laminar  and  turbulent 


diameter  correction  factors  versus  cylinder  diameter  are  presented  In  Fig.  5.2. 


The  difference  In  the  laminar  and  turbulent  curves  stems  from  the  fact  that  the 


laminar  heating  methodology  Is  a function  of  the  cylinder  diameter  to  the  0.5 


power,  whereas  the  turbulent  heating  Is  a function  of  the  diameter  to  0,2  power. 
An  adjustment  Is  made  to  the  corrected  heating  load  obtained  from  Fig.  5.1  for 
the  cylinder  diameter  variation  by  multiplying  this  value  times  the  adjustment 


factor  presented  In  Fig.  5.2.  In  order  to  determine  which  curve,  whether  laminar 


or  turbulent,  is  to  be  used  for  the  adjustment  factor,  use  the  following  procedure. 
For  the  cylindrical  sections  of  both  configuration  types,  use  the  laminar  curve 


for  the  constant  trim  attitude  range  85“  «l  s.  95“.  The  turbulent  curve  Is 
used  for  all  other  trim  attitudes.  For  the  base  section,  use  the  laminar  curve 
for  the  range  175“i,«j  i180“  and  the  turbulent  curve  for  the  remaining  trim 


attitudes. 


In  many  situations  It  may  be  determined  that  the  payload  rolls  during 
reentry.  If  this  Is  the  case  the  heating  load  may  be  reduced  to  the  cylinder 
section  of  the  configuration  by  multiplying  the  factor  Rr/hsA  from  Fig.  5.3 
times  the  diameter  adjusted  heating  obtained  using  Fig.  5.2.  This  factor  Is  to 
be  applied  only  to  cylinder  section  of  the  configuration  since  vehicle  roll  does 


not  reduce  base  region  heating. 


After  the  cylinder  diameter  adjustment  has  been  made,  the  vehicle  surface 


temperature  can  then  be  determined  using  Fig.  5.4.  Fig.  5.4  presents  the 
vehicle  surface  temperature  versus  the  corrected  heating  load.  This  surface 
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temperature  assumes  no  variation  with  temperature  In  the  heat  capacity  of  the 
vehicle's  skin  material.  Although  the  variable  heat  capacity  correction  is  not 
very  significant,  these  corrections  can  be  made  for  the  surface  temperature 
using  Fig.  5.6. 

Although  the  methods  and  procedures  for  determining  surface  temperatures 
on  reentry  vehicles  from  heating  loads  presented  in  this  section  seem  comparatively 
simple,  an  In-depth  analysis  was  performed  to  verify  these  methods.  The  method 
presented  In  the  form  of  Fig.  5.4  is  a "heat  dump  method"  in  which  the  vehicle 
surface,  because  of  the  material  type  and  thickness,  is  assumed  to  react 
instantaneously  to  a heat  load  input  by  an  increase  in  surface  temperature. 

A conduction  heating  analysis  using  a thermal  analyzer  was  performed  to  verify 
the  vallTlity  of  this  assumption.  Detailed  discussions  of  this  analysis  aro 
presented  in  Appendix  D. 

This  section  presents  the  explanations  and  procedures  for  determining 
vehicle  surface  temperatures  from  the  heating  loads  obtained  in  previous 
sections.  The  value  of  the  surface  temperatures  obtained  after  the  variable 
heat  capacity  correction  is  made  Is  the  best  estimate  for  the  surface  temperature 
on  a reentry  vehicle  that  can  be  obtained  using  the  results  presented  in  this 
handbook. 
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Fig.  5.5  Variable  heat  capacity  adjusbuent  to  constant 
heat  capacity  surface  temperature,  T-^. 
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SECTION  6 

PAYLOAD  REENTRY  DYNAMICS 

The  payload  configurations  considered  in  this  analysis  may  either  tumble 
and  oscillate  or  exhibit  normal  stable  flight  at  some  constant  attitude. 

Neither  roll  nor  yaw  type  autorotative  motions  were  coixsidered  in  the  results 
presented  in  this  section  since  the  geometry  and  initial  conditions  of  the  payload 
at  release  did  not  indicate  that  these  types  of  motions  would  occur.  In  Section 
6.1,  ballistic  reentry  (reentry  at  a constant  trim  attitude)  is  discussed  while 
the  oscillating  and  tumbling  trajectories  are  discussed  in  Section  6.2. 

6.1  Ballistic  Reentry  Trajectories 

The  most  Important  concept  behind  ballistic  reentry  is  that  the  tra- 
jectory of  a ballistic  reentry  body  is  a function  only  of  the  body's  ballistic 
coefficient  and  its  initial  altitude.  If  a reentry  body's  weight  to  drag  ratio, 
W/CqA,  defined  as  ballistic  coefficient)  is  known  as  well  as  the 

release  altitude,  then  the  body's  reentry  trajectory  can  be  determined  from 
Figs.  6.1  - 6.9.  These  figures  present  the  reentry  trajectories  for  ballistic 
coefficients  ranging  from  10  to  260  lb/ft*for  the  initial  release  altitudes 
of  75  to  350  statute  miles. 

A given  payload's  geometry  (either  cylinder  or  cone-cylinder)  with  a 
specific  fineness  ratio,  diameter,  and  C.  G.  location  specifies  a constant 
trim  attitude  (See  Section  2)  for  ballistic  reentry.  A number  of  trim  attitudes 
and  ballistic  coefficients  are  presented  in  Table  6,1.  This  information  is  to  be 
used  along  with  the  information  provided  In  Tables  6.2  and  6.3  to  compute  the 
equivalent  ballistic  coefficients  for  the  cases  considered  in  the  analysis. 
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This  coefficient  then  specifies  the  ballistic  reentry  trajectory  for  the  specific 
configuration  type  and  case. 

6.2  Oscillating  and  Tumbling  Trajectories 

In  many  situations,  a certain  C.  G.  location  combined  with  an  Initial 
release  attitude  will  promote  large  and  small  aiiipl Itude  oscll latlons  about  a 
trim  attitude  which  can  be  predicted.  In  many  other  situations  certain  combina- 
tions will  result  In  a complete  tumbling  of  the  vehicle.  Ballistic  trajectories 
can  be  used  in  situations  where  the  vehicle  oscillations  are  small.  Although 
these  are  not  constant  trim  attitude  entries,  the  ballistic  trajectories  closely 
approximate  the  small  amplitude  oscillating  trajectories  without  Inducing 
any  appreciable  error.  In  the  other  situations,  however,  where  the  vehicle 
oscillations  are  large  or  where  the  vehicle  tumbles,  a ballistic  trajectory 
approximation  can  result  In  significant  errors  In  the  actual  trajectory.  Large 
vehicle  oscillations  or  tumbling  result  In  large  variations  In  the  vehicle  drag 
force.  Since  the  descent  velocity  of  a vehicle  Is  dependent  on  Its  drag  force, 
this  results  In  large  variations  In  the  vehicle's  descent  velocl ty  and  thus  Us 
entry  trajectory.  This  section  addresses  the  problem  of  determining  the  reentry 
trajectories  for  the  payload  configurations  considered  In  the  analysis  whose 
variety  of  C.  G.  locations  and  Initial  release  attitudes  result  In  largo  vehicle 
oscillations  and  tumbling  during  reentry. 

A six-degree-of-freedom  trajectory  simulatlcn  computer  code  (Ref.  2) 
was  used  extensively  In  determining  the  effects  of  variations  In  C.  G.  location 
and  initial  release  attitudes  for  the  various  payload  configurations.  This 
computer  code  Is  discussed  in  detail  1n  Appendix  B.  Graphical  presentation 
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Of  the  results  of  the  total  matrix  of  6 DOF  runs  is  not  practical  since  the 
total  number  of  graphs  describing  this  matrix  are  too  numerous.  The  results 
from  a condensed  form  of  the  matrix  of  runs  will  be  presented,  however,  in 
tabular  form  in  Tables  6.2  and  6.3.  From  the  6 DOF  trajectory  runs,  an  equivalent 
ballistic  coefficient  was  obtained  in  situations  witere  the  trajectory  varies 
from  a constant  trim  entry  ballistic  trajectory.  The  percent  of  variation  of 
equivalent  ballistic  coefficients  from  the  constant  trim  ballistic  coefficients 
were  determined.  These  results  are  shown  in  Table  6.2.  While  Table  6.1  gives 
the  ballistic  coefficients  for  normal  stable  constant  trimmed  entries  of  the 
payload  configurations  and  cases  being  analyzed,  Table  6.2  gives  the  variations 
from  these  coefficients  when  initial  vehicle  attitudes  combined  with  various 
C.  G.  offsets  are  considered.  If  this  variation  is  positive,  zero,  or  negative, 
this  Indicates  that  tho  equivalent  ballistic  coefficient  for  the  specific  case 
is  either  greater  than,  the  same,  or  less  than  the  constant  trim  ballistic 
coefficient  by  the  magnitude  of  the  percent  shown.  Table  6.3  specifies  the 
payload  dynamics  as  to  whether  it  tumbles,  oscillates,  or  is  trimmed  during  entry. 
From  Table  6.3,  it  is  possible  to  determine  if  the  vehicle  is  reentering 
forward  or  aft  end  first.  If  a payload  is  oscillating,  it  is  also  possible  to 
determine  the  degree  of  severity  of  these  oscillations,  i.e.,  whether  the  payload 
oscillations  have  small  amplitudes  or  whether  the  amplitudes  are  large  and  unstable. 

Perhaps  the  main  criticism  in  presenting  data  in  a tabular  fashion  as  in 
Tables  6.2  and  6.3  is  how  the  data  is  supposed  to  be  applied  when  the  payload 
geometry,  C.  G.  location,  or  initial  release  attitudes  vary  from  the  cases 
presented.  It  is  possible,  for  instance,  although  unlikely,  that  a cone-cylinder 
configuration  of  fineness  ratio  4 1/2  may  exhibit  reentry  dynamics  which  are 
drastically  different  from  the  dynamics  exhibited  by  the  cone-cylinder  geomotries 
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with  fineness  ratios  of  4 and  5 which  are  presented  in  the  tables.  Of  course, 
this  will  always  be  a possibility  since  there  is  no  way  to  increase  the  6 DOF  run 
niatrix  to  include  all  fineness  ratios.  It  is  possible  to  observe  certain  trends 
in  Tables  6,2.  and  6.3  when  considering  different  fineness  ratios.  Notice  that, 
in  general,  C.  G.  locations  aft  of  the  volume  C.  G.  promote  reduced  ballistic 
coefficients  for  the  cone-cylinder  configurations  whereas  forv/ard  C.  G.  locations 
tend  to  Increase  the  ballistic  coefficients.  This  is,  of  course,  with  the 
exception  of  the  i/d  equal  to  3 configuration  which  is  essentially  identical 
to  stable  ballistic  reentry  no  matter  what  combination  of  C.  G.  offset  or  initial 
release  attitude  is  being  considered.  With  the  exception  of  the  1055  forward 
C.  G.  offset  and  the  180"  initial  release  attitude  combination,  C.  G.  locations 
from  the  volumetric  C.  G.  to  10%  forward  for  the  initial  release  attitudes  of 
150"  to  180"  in  general,  promote  the  largest  increase  in  ballistic  coefficients 
for  the  cone-cylinder  geometries.  The  reason  for  the  large  increase  in  ballistic 
coefficients  for  these  combinations  is  due  primarily  to  the  large  increase  in 
descent  velocity  which  results  from  large  decreases  in  drag  force.  Large  variations 
in  the  vehicle  attitude  during  a vehicle  oscillating  entry  allow  the  vehicle  to 
fly  at  some  attitudes  where  the  drag  force  is  very  small,  thus,  creating  large 
descent  velocity. 

Similarly,  decreases  in  ballistic  coefficients  can  be  understood  when 
considered  in  this  manner.  A number  of  other  trends  can  be  established  and 
rationalized  for  both  the  cylinder  and  cone-cylinder  configurations.  The  data, 
although  not  necessarily  linear  between  fineness  ratios,  C.  G.  offsets,  and 
initial  release  attitudes,  is  consistent  and  can  be  applied  in  situations 
when  a flight  case  and  a case  presented  in  Table  6.2  do  not  coincide.  When  this 
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situiitlon  occurs,  first  dotoriuine  tiu)  bounds  for  tho  fli.jitt,  cnso.  These  include 
the  confitjuration  typo,  and  tho  upper  and  lower  bounds  on  fineness  ratio, 

C.  G.  location,  and  release  attitude.  Tho  percertt  of  variation  of  the  ballistic 
coefficient  is  then  determined  for  these  bounds,  Tiie  percent  of  variation 
should  be  weighted  more  strongly  to  tho  results  which  are  closest  to  the  fliglit 
case.  This  weighting  can  be  performed  by  using  linear  interpolation  or  by 
considering  the  case  in  more  detail  and  determining  a weigl^ting  factor  based 
on  engineering  judgement. 

The  equivalent  ballistic  coefficients  discussed  in  this  section  relate  only 
to  tho  dynamics  of  tho  vehicle.  They  are  not  to  lie  coni'used  with  or  applied 
to  tha  heating  environment  for  a vehicle  during  an  oscillating  or  tumbling 

reentry.  In  addition  to  ti>e  figures  and  tables  presented  in  this  section  to 

«• 

determine  trajectories  for  reentry  vehicles,  rig.  6.10  has  been  provided  in 
order  to  assist  the  engitKior  in  the  aerial  recovery  .of  payloads.  The  maximum 
decent  velocity  of  a vehicle  is  presented  "orsus  ballistic  coefficient  and  release 
altitude  for  two  recovery  altitudes.  This  should  provide  tl\e  engineer  with  a 
quick  indicator  of  possible  dynamic  problems  that  could  arise  and  effect  the 
siKcessful  recovery  of  Lite  vehicle. 
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Geometry  of  Cone- Cylinder  Confiyuratinns 
|— £ ^ 

~ I ^ Cody  1 , Jl/d  =>  7 


C 

C 


e 


Body  2,  £/d  ■ S 
Body  3,  i/d  = 4 
j Body  4,  i/d  ^ 3 

6 Volumetric  C.  G.  Location 


Conditions:  Density  assumed  constant  and  ® 23.23  Ib/ft*;  diameter  » 2.5  ft. 


6o(1y 

C.G. 

Location 

Trim 

Attitude  (aj)* 

Weight  (lb) 

Drag 

Coef.  (C[j) 

Ballistic 
Coef.  (Ib/ft*) 

1 

lOX  F 

37 

1614 

2.35 

140 

2 

m F 

44 

1044 

2.40 

89 

3 

m F 

45 

759 

1.95 

79 

4 

103!  F 

30 

474 

0.7 

138 

1 

'53!  F 

64 

1614 

7.2 

46 

2 

S%  F 

73 

1044 

5.7 

37 

3 

S%  F 

85 

759 

5.0 

31 

4 

53!  F 

180 

474 

1.6 

60 

1 

2.63!  F 

90 

1614 

10.1 

33 

2 

2.53!  F 

107 

1044 

6.3 

34 

3 

2.53!  F 

147 

759 

2.6 

59' 

4 

2.5%  F 

180 

474 

1.6 

60 

1 

Vol  C.G. 

115 

1614 

8.3 

40 

2 

Vol  C.G. 

137 

1044 

3.7 

57 

3 

Vol  C.G. 

180 

759 

1.6 

97 

4 

Vol  C.G. 

180 

474 

1.6 

60 

1 

2.53!  Aft 

138 

1614 

4.6 

71 

2 

2.53!  Aft 

158 

1044 

2.1 

101 

3 

2.53!  Aft 

180 

759 

1.6 

97 

4 

2.53!  Aft 

180 

474 

1.6 

60 

1 

531  Aft 

156 

1614 

2.5 

132 

2 

53!  Aft 

168 

1044 

1.8 

118 

3 

531  Aft 

180 

759 

1.6 

97 

4 

0%  Aft 

180 

474 

1.6 

60 

1 

103!  Aft 

168 

1614 

1.7 

193 

2 

10%  Aft 

18' 

1044 

1.6 

133 

3 

1031  Aft 

loO 

759 

1.6 

97 

4 

103!  Aft 

180 

474 

1.6 

60 

Table  6.1a  Information  for  determining  the  Ballistic  Trajectories 
for  particular  cone-cylinder  configurations 


GEOMETRY  OF  CYLINDER  CONFIGURATIONS 


1_. 

Body  5,  £/d  =2.5 
Body  6,  Vd  ■=  4 


Body  7,  Vd  = 5 


e Volumetric  C.  6.  Location 


Conditions:  Density  assumed  constant  and  = 23.23  Ib/ft  ; diameter  = 2.5  ft. 


Body 

C.G. 

Location 

« 

Trim 

Attitude  (ut)® 

Weight  (lb) 

Drag 

Coef.  (Cp) 

Ballistic 

Coef. 

5 

10?i  F 

59 

499 

3.6 

28 

6 

1051  F 

48 

980 

4.0 

50 

7 

10%  F 

42 

1225 

3.9 

64 

5 

5%  F 

74 

499 

4.2 

24 

6 

551  F 

66 

930 

5.9 

34 

7 

5%  F 

62 

1225 

6.6 

38 

5 

2.5%  F 

82 

499 

4.3 

24 

6 

2.5%  F 

70 

980 

6.5 

31 

7 

2.5%  F 

75 

1225 

8.0 

31 

5 

Vol  C.G. 

90 

499 

4.4 

23 

6 

Vol  C.G. 

90 

980 

6.7 

30 

7 

Vol  C.G. 

90 

1225 

8.4 

30 

5 

2.5%  Aft 

98 

499 

4.3 

24 

6 

2.5%  Aft 

102 

980 

6.5 

31 

7 

2.5%  Aft 

105 

1225 

8.0 

31 

5 

5%  Aft 

106 

499 

4.2 

24 

6 

5%  Aft 

114 

980 

5.9 

34 

7 

5%  Aft 

118 

1225 

6,6 

38 

5 

10%  Aft 

121 

499 

3.6 

28 

6 

10%  Aft 

132 

980 

4.0 

50 

7 

10%  Aft 

138 

1225 

3.9 

64 

Table  6.1b 


Information  for  determining)  the  ballistic  trajectories  for  particular 
cylinder  configurations. 
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Fig.  6.1  Ballistic  trajectories  for  an  initial  release  altitude  (H^-)  of  75  miles. 
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VELOCITY,  FT/SEC  TIME  IN  SECONDS 

Fig.  6.6  Ballistic  trajectories  for  an  initial  release  altitude  (H.j)  of  225  miles. 
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Fig.  6.7  BalTistfc  trajectories  for  an  initial  release  altitude  (H^ ) of  Z50  miles. 
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6.9  Ballistic  trajectories  for  an  initial  release  altitude  (H.)  of  350  miles. 
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SECTION  7 

CONVECTIVE  HEATING  RATES 

Convection  heating  rates  versus  altitude  and  velocity  for  angles  of 
attack  ranging  from  90  to  180  degrees  are  presented  In  Figs.  7>1  - 7.10. 

These  heating  rates  are  applicable  to  the  base  and  cylindrical  sections  of 
both  the  cylinder  and  cone-cylinder  configurations.  A method  for  using  the 
plots  for  the  cylinder  heating  rates  to  compute  heating  rates  for  the  cone 
of  the  cone-cylinder  configurations  Is  given  In  Table  7.1.  Detailed  dis- 
cussions of  the  convection  heating  methodologies  used  In  the  analysis  are 
presented  1n  Appendix  C. 

These  heating  rates  are  presented  In  this  section  In  order  to  familiarize 
the  design  engineer  with  heating  rate  environments  for  reentry  payloads.  It  Is 
entirely  feasible  that  certain  pieces  of  hardware  on  board  a reentry  payload 
could  be  heating  rate  sensitive,  1.e.,  the  piece  of  hardware  Is  unable  to  withstand 
an  extremely  high  heating  rate  for  a short  time  period.  Certain  types  of  In- 
strumentation could  possibly  fall  Into  the  category  of  being  rate  sensitive 
when  subjected  to  high  reentry  heating  rates. 

There  Is  an  additional  benefit  to  the  design  engineer  In  being  able  to 
determine  the  heating  rate  history  for  a vehicle.  If  the  altitude,  velocity, 
time  and  angle  of  attack  history  are  known  for  a reentry  payload,  an  Integrated 
heating  load  calculation  can  be  performed.  A surface  temperature  for  the  reentry 
payload  can  then  be  determined  from  the  Integrated  heating  load.  If  the  vehicle 
Is  reentering  at  a constant  angle  of  attack,  a heating  rate  at  specific  time 
points  con  be  obtained  by  superimposing  the  payload  trajectory  onto  the  correct 
angle  of  attack  figure  (Figs.  7.1  - 7.10).  A nunu.,  cal  Integration  of  the  heating 
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rate  over  a certain  time  Interval  can  be  performed  1n  order  to  determine  an 
Integrated  heating  load.  This  is  essentially  the  same  method  by  which  heating 
loads  were  comouted  in  Section  3 of  this  handbook  except,  of  course,  the 
methods  were  computerized  in  order  to  cut  down  on  the  amount  of  work  required 
and  Increase  the  accuracy  of  the  calculation.  If  the  payload  is  not  trimmed 
during  entry,  but  rather  oscillating  or  tumbling,  an  integrated  load  can  still 
be  determined  in  similar  fashion  provided  the  vehicle  attitude  is  known  at  a 
number  of  selected  tinw  points  during  the  peak  heating  portion  of  the  entry 
trajectory. 

The  heating  rate  history  of  a reentry  vehicle  can  also  be  useful  in  other 
areas.  The  conduction  heating  from  the  surface  to  various  structural  members 
inside  the  vehicle  reduces  the  surface  temperature.  This  reduction  in  surface 
temperature  can  be  determined  using  a thermal  analyzer  conduction  program 
once  the  heating  rate  history  and  a model  of  the  structural  members  are  obtained. 
This  type  of  analysis  is  recommended  when  the  handbook  results  indicate  excessive 
surface?  temperatures  beyond  the  safe  limits  of  the  material.  The  analysis  is 
difficult,  time  consuming,  and  requires  a thorough  knowledge  of  conduction 
heat  transfer  and  thermal  analyzer  computer  programs.  Even  though  this  is  the 
case,  the  analysis  should  be  considered  in  these  situations., 

During  a vehicle  reentry  there  is  an  altitude  at  which  the  aerodynamic 
heating  ends  and  cooling  begins.  Fig.  7.11  presents  this  altitude  versus 
ballistic  coefficient  for  three  release  altitudes  of  the  vehicle.  The  altitudes 
presented  in  Fig.  7.11  provide  cutoff  limits  for  the  heating  rate  analyses. 

These  same  limits  were  used  in  determining  the  loads  presented  1n  Section  3. 

The  heating  rate  plots  (Figs.  7.1  to  7.10)  are  based  on  a 24-1nch  reference 
diameter.  The  flow  regime  transition  from  laminar  to  turbulent  flow  is  built 
ill  the  filots  since  the  transition  criterion  is  independent  of  diameter.  The 
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heating  rate  from  the  figures  (Fig.  7.1  to  7.10)  may  be  corrected  for  diameter 
changes  using 

*1  ’ ‘•Fig  ^24  inches  ) 

where 

qf^g  « heating  rate  from  Figs.  7.1  - 7.10 

a ■ 0.5  laminar 
a ■ 0.2  turbulent 

In  order  to  determine  a value  for  the  constant,  a,  In  the  above  expression 
the  flow  regime  must  be  established.  The  free  stream  Reynolds  number  per  foot 
(Re«yft)  will  be  used  to  Indicate  whether  the  flow  Is  laminar  or  turbulent  at 
the  particular  altitude  velocity  point  where  the  heating  rate  Is  desired.  The 
expression  for  Re*/ft  1$ 

Re«/ft  ■ 

where 

p«o  ■ free  stream  density  (1bm/ft*) 

Vo#  ■ free  stream  velocity  (ft/sec) 

Poo  ■ free  stream  viscosity  (Ibm/ft-sec) 

The  free  stream  velocity  In  the  Reynolds  number  expression  Is  known.  The 
free  stream  density  can  be  obtained  at  the  given  altitude  from  Fig.  7,15.  Ir 
order  to  determine  m«,  Sutherland's  law  for  viscosity  will  be  used  along  with 
the  free  stream  temperature  from  Fig.  7.14.  The  expression  for  Sutherland's 
law  Is 

Po#  • 7.3036  X 10“’  r«o‘“Vw:g'  ^b"’/ft-sec  (T«,  In  ®R) 

Once  Re«/ft  has  been  computed  the  flow  regime  (laminar  or  turbulent)  may  be 
determined  from  Fig.  7.12.  In  general#  the  flow  regime  will  be  laminar  for 
altitudes  greater  than  300K  feet  and  turbulent  for  altitudes  less  than  50K  feet. 
Certain  angles  of  attack  will  promote  laminar  or  turbulent  flow  regardless  of 
the  altitude  and  velocity'. 
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Cylinder;  85  < < 95  Always  laminar 

Base:  c»i-  > 175  Always  laminar 

If  angle  of  attack,  Oj,  Is  less  than  90  degrees,  the  following  procedure 
Is  used 

Cylinder;  If  oj  <90,  compute  the  heating  using  180  - aj 
Base:  If  a-p  < 90,  compute  the  heating  using  ay  •>  90. 

Once  the  heating  rate  Is  computed  using  the  preceding  procedure,  a surface 
temperature  correction  may  be  required.  The  results  presented  1n  Fig.  7.1  to 
Fig.  7.10  are  for  a wall  temperature  of  660**R.  The  heating  rate  may  be 
corrected  for  wall  temperature  by  using  the  following  equation. 

where 

• T-  ♦ 5?;^  ex) 

2CpgJ  - 12019.2  (ft*/sec**’R) 

R ■ Recovery  factor  from  Fig.  7.13 

T«  ■ Free  stream  temperature  from  Fig.  7.14 

Ts  « Surface  temperature  ("R) 

Thus,  Fig.  7.1  to  7.10  may  be  used  by  making  angles  of  attack,  diameter,  and 
wall  temperature  adjustments  to  the  reference  conditions  for  a large  variety  of 
problems. 
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TABLE  7.1 

CONE  HEATING  RATES 


The  procedure  for  computing  heating  rates  to  the  cone  of  a cone- 
cylinder  configuration  is  as  follows: 

(1)  Determine  the  angle  of  attack  of  the  cone's  windward  surface 

cone  cyl  c 

“cyl  " attack  of  cylinder 

Oj.  ■ Semi -vertex  angle  of  the  cone 

NOTE:  If  > ISO*,  ■ 180' 
cone  cone 

(2)  Determine  the  heating  rata  for  a cylinder  at  from  F1g.  7.1 
to  7.10. 

(3)  Determine  the  cone  radius  at  the  mldlcngth  of  the  cone»  and  cone- 
cylinder  radius  ratio.  The  reference  radius  for  the  cylinder  Is 
'"cyl  " Inches  In  Fig.  7.1  to  7.10. 


(4) 


Correct  the  heating  rate  for  the  rad'lus  effect. 


where  a ■ .5  laminar 
a ■ .2  turbulent, 

and  where  the  type  of  flow  can  be  determined  from  Fig.  7.12  based 
on  altitude  and  velocity.  If  the  angle  of  attack,  Is  between  65 
and  95  degrees  the  flow  Is  laminar  regardless  of  the  altitude  and 
velocity. 
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Ffg.  7.6  Convection  heating  rates  for  aj  = 140  degrees 
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Fig.  7.7  Convection  heating  rates  for  aj  = 150  degrees. 
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Fig.  7.8  Convection  heating  rates  for  ay  = 160  degrees 
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End  of  Positive  Heat  (T  =660“^ 
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F1g.  7.12  Flow  vmjline  transition  criterion  tor  nono-cyl  Indors  at 
largo  angle  of  attack. 
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SECTION  8 
REENTRY  GUIDELINES 

In  general,  reentry  problems  are  coupled  problems  In  that  It  is  Impossible 
to  separate  and  analyze  only  the  dynamics  or  heating  aspect  without  considering 
the  other.  This  Is  especially  true  when  considering  missions  with  high  Initial 
release  altitudes.  The  higher  the  release  altitude  the  more  coupled  the  dynamics 
and  heating  becomes.  This  Is  the  reason  that  certain  restraints  and  guidelines 
must  be  followed  when  considering  a change  In  the  design  of  a reentry  vehicle  to 
solve  a dynamics  or  heating  problem. 

8.1  Mission  and  Configuration  Dependent  Effects 

There  are  a number  of  mission  and  configuration  dependent  effects  associated 
with  the  dynamics  and  heating  of  a reentry  payload.  These  Include  the  effects 
of  (1)  C.  G.  location,  (2)  release  attitude,  (3)  release  altitude,  (4)  geometric 
configuration,  (5)  skin  thickness,  (6)  weight,  and  (7)  external  protuberances. 

The  mission  usually  dictates  which  effects  can  and  cannot  be  controlled.  This 
In  turn  leaves  the  engineer  In  many  cases  with  very  little  flexibility  In 
designing  around  a potential  dynamics  or  heating  problem.  The  Initial  release 
altitude,  configuration  type,  weight,  and  certain  external  protuberances  are 
generally  felt  to  be  mission  dependent.  The  release  attitude  and  the  C.  G. 
location  of  the  vehicle  are  considered  to  be  more  flexible  and  not  totally  mission 
dependent.  If  the  weight  requirements  of  the  mission  are  not  too  stringent, 
there  could  possibly  be  additional  flexibility  from  the  design  standpoint  In 
being  able  to  Increase  the  skin  thickness  of  the  vehicle  In  certain  areas.  In 
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many  situations,  a potential  dynamics  or  heating  problem  can  be  relieved  by 
modifying  the  C.  G.  location,  controlling  the  vehicle  release  attitude,  or 
varying  the  vehicle  skin  thickness  provided  certain  guidelines  are  followed  when 
doing  so. 

As  previously  discussed,  a reentry  problem  Is  a coupled  problem  In  that  any 
change  In  the  reentry  dynamics  of  the  vehicle  can  and  usually  does  cause  a 
corresponding  change  In  the  heating  environment.  The  converse  Is  also  true. 
Because  of  this,  the  dynamics  and  heating  problems  should  be  solved  simultaneously 
by  an  Iterative  type  procedure.  In  order  to  demonstrate  this  procedure,  consider 
the  case  In  which  a dynamics  problem  has  been  Indicated  by  the  results  In  the 
previous  sections.  A shift  In  the  C.  G.  location  Is  a possible  solution  to 
problems  since  each  vehicle  has  a C.  G.  location  which  promotes  a vehicle  trim 
attitude  at  which  maximum  drag  occurs.  This  results  In  a minimum  descent  velocity 
for  the  vehicle.  If  this  velocity  Is  acceptable  for  the  recovery  system  being 
used  then  the  heating  environment  for  the  vehicle  should  be  determined  based  on 
the  maximum  drag  attitude.  If  the  resultant  heating  environment  does  not  create 
any  difficulties,  then  the  total  reentry  problem  has  been  solved  by  modifying 
the  vehicle  C.  G.  location.  In  the  event  that  a severe  heating  environment 
results,  then  Iterate  until  an  optimum  C.  G.  location  Is  determined  which  will 
satisfactorily  solve  both  problems.  For  a number  of  missions,  particularly  the 
higher  altitude  missions,  there  will  be  no  optimum  solution  to  the  combined 
problem  of  reentry  dynamics  and  heating.  When  a C.  G.  location  Is  varied  to  solve 
a dynamics  problem,  a severe  heating  problem  will  be  created.  Conversely, 
proposed  solutions  to  the  heating  problem  will  create  severe  dynamic  problems. 

When  this  situation  arises,  serious  consideration  should  be  given  to , performing 
a complete  theoretical  analysis  for  the  particular  mission  including  a 6 DOF 
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trajectory  analysis  and  a tumbling  heating  analysis.  These  results  will  be 
more  accurate  than  the  generalized  results  obtained  from  the  handbook.  If  these 
results  still  indicate  severe  dynamic  and  heating  problems,  then  another  option 
should  be  considered  such  as  inducing  some  effect  by  Initiating  a roll,  yaw,  or 
tumble  motion  to  the  vehicle.  The  induced  effects  will  be  discussed  in  the 
following  section.  If  it  is  not  possible  to  induce  certain  effects,  or  if  it  is 
generally  ffelt  that  these  effects  will  not  satisfactorily  solve  the  problem, 
then  serious  consideration  should  be  given  to  changing  the  vehicle  design  or 
the  method  used  in  recovering  the  vehicle. 

The  surface  temperature  of  a reentry  payload  can  be  reduced  by  increasing  the 
skin  thickness.  If  the  skin  thickness  is  increased  in  a uniform  manner,  or  such 
that  the  vehicle  C.  G.  location  is  not  changed,  this  will  cause  an  increase  in 
the  descent  velocity  of  the  vehicle  due  to  the  added  weight.  The  attitude  dynamics 
of  the  payload  will  remain  essentially  unchanged.  in  general,  increases  in 
the  skin  thickness  of  a payload  such  that  the  vehicle  C.  6.  location  remains 
unchanged  will  have  the  net  effects  of  decreasing  the  surface  temperature  and 
increasing  the  maximum  descent  velocity  of  the  vehicle.  If  the  vehicle  skin 
thickness  is  increased  in  a fashion  such  that  it  causes  a shift  in  the  C.  G. 
location,  this  changes  the  vehicle  attitude,  descent  velocity  and  also  the  heating 
environment.  When  an  Increase  in  the  vehicle  skin  thickness  is  being  considered, 
a complete  analyses  should  be  performed  using  the  results  presented  in  the 
handbook  in  order  to  determine  the  net  effect  of  this  Increase  from  both  a 
dynamics  and  heating  standpoint. 

The  results  presented  in  this  handbook  apply  specifically  to  clean  skin 
reentry  vehicles  (i.e.  vehicles  with  no  external  protuberances).  Even  though  this 
is  the  case,  a certain  amount  of  Information  concerning  protuberances  will  be 
provided  in  this  section  since  the  presence  of  external  protuberances  on  a 
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reentry  vehicle  can  drastically  alter  Its  reentry  characteristics.  Protub- 
erances are  pluced  on  payloads  for  many  functions.  Some  are  required  In  order 
for  the  payload  to  perform  Its  mission  as  would  be  the  case  of  external  equip- 
ment used  In  data  acquisition.  Other  types  of  protuberances  are.  placed  on  the 
payload  to  Induce  some  desired  dynamic  effect  of  the  vehicle.  These  types  of 
protuberances  will  be  discussed  In  the  following  section.  When  external  pro- 
tuberances are  attached  to  a payload,  the  aerodynamic  characteristics  of  that 
payload  can  be  significantly  changed.  The  extent  to  which  these  characteristics 
are  changed  depends  upon  the  size,  weight,  and  location  of  these  protuberances. 

In  addition  to  changing  the  aerodynamic  characteristics,  they  can  also  create 
additional  heating  problems  because  of  the  Increased  localized  heating  In  areas 
on  and  around  the  protuberance.  If  external  protuberances  are  required  for 
a payload,  a symmetrical  circumferential  distribution  of  these  protuberances 
Is  desired  If  at  all  possible.  If  the  protuberances  are  placed  on  the  payload 
In  an  unsymmetric  fashion,  the  vehicle  will  tend  to  reach  an  equilibrium 
position  such  that  the  protuberances  are  on  the  leeward  side  of  the  vehicle 
away  from  the  flow.  This  Is  a disadvantage  in  many  circumstances.  For  Instance, 

If  an  Initial  roll  rate  was  Imposed  on  the  payload  at  separation  In  order  to 
relieve  the  heating  load  during  reentry,  an  unsymiiietrlcal  protuberance  distri- 
bution could  possibly  damp  out  the  vehicle  roll.  This  would  cause  one  side  of 
the  vehicle  to  be  subjected  to  maximum  stagnation  line  heating  Instead  of 
distributing  the  load  which  was  the  original  aim  In  trying  to  Induce  vehicle 
roll.  An  unsymmetrical  distribution  of  protuberances  can  also  result  In 
extremely  complex  vehicle  dynamics.  Various  combinations  of  roll,  yaw  and  tumble 
motions  could  result  which  would  have  significant  effects  on  the  heating  envi- 
ronment and  the  eventual  recovery  of  the  vehicle.  External  protuberances  can  create 
a whole  host  of  other  problems  In  the  recovery  phases  of  a reentry  payload. 
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In  general,  the  effects  of  protuberances  with  heights  above  the  external  surface 
less  than  approximately  0.5<  of  the  vehicle  diameter  can  be  neglected.  The 
effects  of  external  protuberances  larger  than  this  should  be  determined  using 
a separate  analysis  for  the  particular  configuration  and  mission. 

8.2  Induced  Effects  for  Reentry  Payloads 

Certain  induced  effects  can  be  applied  to  control  the  dynamics  of  a 
vehicle  in  order  to  obtain  a satisfactory  solution  to  a heating  or  dynamics 
problem.  Roll,  yaw  and  tumble  motions  can  be  Imparted  to  the  payload  at 
separation.  These  types  of  nK)tions  can  also  result  from  the  addition  of  selected 
protuberances  placed  in  specific  locations  on  the  vehicle.  Roll,  yaw  and  tumble 
motions  and  the  methods  used  to  initiate  and  sustain  these  types  of  motions  are 
the  subjects  of  this  section. 

If  a payload  can  sustain  a high  roll  rate  during  reentry,  the  heating  load 
to  the  cone  and  cylinder  sections  of  the  vehicle  types  considered  can  be  signifi- 
cantly reduced.  This  roll  rate  should  be  a minimum  of  1/2  revolution  per 
second  as  initial  roll  rate  can  be  imparted  to  the  payload  at  release.  This 
roll  rate  should  continue  and  the  heat  load  reduced  provided  this  motion  is 
not  damped  out  by  external  protuberances.  Roll  rates  can  also  be  initiated  and 
sustained  during  reentry  by  attaching  aerodynamic  vanes  or  rotors  to  the  external 
surface  of  the  reentry  payload.  Reference  l presents  data  taken  in  the 
University  of  Maryland  wind  tunnel  which  Indicates  that  high  roll  rates  can  be 
obtained  by  attaching  low  aspect  ratio  fins  or  rib-type  rotors  to  the  external 
surface  of  the  vehicle.  No  matter  what  method  is  selected  to  initiate  and 
sustain  a roll  rate,  rolling  a reentry  vehicle  is  an  excellent  way  to  reduce 
the  heating  load  on  the  cylinder  and  cone  sections  of  the  configurations  being 
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considered.  The  heating  analyses  presented  In  thl:*  handbook  assumes  that  no 
roll  iDotlon  has  been  applied  to  the  vehicle.  Accordingly,  the  heating  load  Is 
calculated  for  the  stagnation  line.  Since  the  precise  roll  orientation  of  the 
body  Is  not  known  during  reentry,  the  entire  vehicle  must  be  designed  for  this 
maximum  condition.  If,  however,  a roll  rate  can  be  Initiated  and  sustained 
during  reentry,  the  heating  load  may  be  reduced  by  considering  the  Information 
provided  In  Fig,  6.3.  This  figure  gives  the  roll  averaged  heat  transfer 
coefficient  ratloed  to  the  stagnation  line  value.  This  ratio  can  be  used  as 
a multiplication  factor  times  the  stagnation  line  load  computed  from  the  previous 
sections.  It  should  be  noted  that  rolling  the  vehicle  will  not  yield  a heating 
relief  to  the  base.  It  Is  more  desirable  to  Induce  a rolling  motion  as  opposed 
to  a yaw  or  tumble  motion  In  that  this  typo  of  motion  Is  easier  to  Induce  and 
sustain.  It  should  be  noted  also  that  a roll  motion  will  have  no  effect  on 
the  descent  velocity  of  the  vehicle  and  should  be  considered  only  to  reduce 
the  heating  environment. 

An  excellent  way  to  reduce  the  descent  velocity  of  a reentry  vehicle  is 
by  initiating  a yawing  or  flat  spinning  autorotative  motion  of  the  vehicle. 

This  motion  will  allow  tlie  vehicle  to  maintlan  a high  angle  of  attack  during 
reentry.  This  high  angle  of  attack  will  create  o maximum  drag  force  which  will 
significantly  reduce  the  descent  velocity  and  old  In  the  recovery  of  the  vehicle. 
Recently,  experiments  have  shown  that  flaps  and  fins  strategically  placed  on 
cylinder  and  cone-cyl1nder  type  configurations  allow  a vehicle  to  sustain 
this  type  of  motion  during  reentry.  Even  though  these  types  of  motions  could 
greatly  aid  In  the  recovery  of  a payload,  consideration  must  be  given  to  the 
heating  environment  which  results  from  the  vehicle  flying  at  these  angles  of 
attack  during  reentry.  Consideration  must  also  be  given  to  the  Increase  In 
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localized  heating  due  to  addition  of  flaps  and  fins  to  the  external  surface 
of  the  vehicle. 

Another  method  of  reducing  the  heating  environment  Is  to  Initiate  a 
tumbling  motion  of  the  vehicle.  If  no  vehicle  roll  Is  considered,  a tumbling 
motion  results  In  distributing  the  heating  load  to  both  sides  of  the  vehicle 
Instead  of  only  heating  one  side  which  Is  the  case  for  a constant  trim  entry 
with  no  roll.  The  descent  velocity  of  a vehicle  cannot,  In  general,  be  reduced 
significantly  with  tumble  unless  the  trim  attitude  of  the  vehicle  Is  such  that 
minimum  or  near  minimum  drag  force  results.  A tumble  motion  can  create  addi  • 
tional  problems  from  the  standpoint  of  parachute  deployment  for  recovery. 

These  possibilities  must  be  considered  before  Initiating  a tumble  motion  to 
relieve  a heating  problem. 

8.3  Material  Selection 

The  handbook  procedures  and  tables  for  computing  a surface  temperature  In 
Section  7 are  based  on  using  aluminum  as  the  vehicle  skin.  If  another  material 
Is  used  or  the  heating  load  Is  sufficiently  high  to  require  a different  material, 
an  alternate  procedure  Is  required.  In  these  events,  the  surface  temperature 
may  be  estimated  using  the  following  equation. 

Tj  - AT  + 660  8.1 

AT  » q/pC6  8.2 

Where 

Q ■ heating  load  C ■ specific  heat 

p ■ material  density  6 ■ material  thickness 

Equations  8.1  and  8.2  are  applicable  only  to  homogeneous  materials  and  should 
not  be  applied  to  any  type  of  nonhomogeneous  material  such  as  a honeycomb 
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structure  for  Instance.  This  relation  also  assumes  that  the  material  Is 
sufficiently  thin  to  respond  Instantaneously  to  the  heating  rate  Input  and 
negligible  surface  reradlation. 


TABLE  8.1 

TYPICAL  HEAT  SINK  MATERIALS 

Material 

Melting 

PoInt’R 

Density 

Ibm/ft 

Specific 

Heat 

Btu/lbm'’R 

(Tw)max 

•R 

A1 umlnum 

1680.1 

168.7 

0.215 

1020 

St.  Steel 
(301) 

Magnesium 

3030 

494.2 

0.108 

2310 

1662 

108.8 

0.25 

1002 

Copper 

2441.1 

559. 

0.092 

1781 

Alternate  heat  sink  materials  are  listed  In  Table  8.1  If  the  material 
thickness  for  any  of  those  listed  In  Table  8.1  exceeds  0.50  Inches • the  temper- 
ature lag  through  the  surface  may  be  significant.  This  would  violate  the 
assumptions  Inherent  In  using  Equations  8.1  and  8.2.  A transient  analysis  Is 
then  required. 

8.4  Initial  Vehicle  Surface  Temperature 

The  maximum  surface  temperature  of  a reentry  vehicle  Is  dependent  upon  Its 
Initial  surface  temperature  due  to  ascent.  A heating  analysis  was  performed 
on  an  Arles  I payload  (cone-cylinder  configuration)  for  three  typical  ascent 
trajectories.  These  trajectories,  obtained  from  Ref.  3,  were  typical  of 
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Heating 

T,/3CR) 

T,^4CI» 

Ti/jCR) 

Mission 

Region 

Light 

Cylinder 

723. 

637. 

590. 

Intermediate 

Cylinder 

710. 

630. 

586. 

Cone 

863. 

723. 

635. 

Hemisphere 

1002. 

790. 

665. 

Heavy 

Cylinder 

696. 

623. 

583. 

Table  8.2  Initial  surface  temperatures  on  an  Arles  I payload  for  three 
ascent  trajectories 

The  three  typical  trajectories  are  designated  as  light,  heavy,  and  Inter- 
mediate missions.  The  light  mission  considers  higher  ascent  velocities  and 
heating  while  the  heavy  mission  considers  ascent  velocities  and  heating  which 
are  significantly  lower.  The  Intermediate  trajectory  results  are  between 

these  two  extremes.  The  Initial  surface  temperatures  denoted  at  Two*  T, ... 

l/o  i/4 

and  Tyg  ^0'"  wall  thicknesses  of  1/8,  1/4  and  1/2  Inches  respectively. 

It  is  apparent  from  Table  8.2  that  the  Initial  surface  temperature  Is  more 
highly  dependent  on  the  vehicle  skin  thickness  than  Its  ascent  trajectory. 

From  Table  8.2  It  Is  also  apparent  that  the  Intermediate  results  provide 
adequate  nominal  estimates  of  the  Initial  surface  temperatures  attained  during 
ascent  for  each  section  of  the  reentry  vehicle. 

The  surface  temperatures  provided  In  Table  8.2  consider  no  cooling  due 
to  radiation  losses  from  the  vehicle  surface  that  would  occur  when  the  vehicle 


H 

{ 

i: 
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Is  outside  the  sensible  atmosphere.  Calculations  were  performed  which  Indicate 
that  these  losses  are  negligible.  The  ascent  analysis  also  considers  no 
reduction  In  the  Initial  surface  temperature  due  to  conduction  heat  transfer 
from  the  surface  to  various  structural  members  and  bulk  heads  located  Inside 
the  vehicle.  The  reduction  in  the  initial  as  well  as  the  niaxiinuiii  surface 
temperature  can  be  considerable  dependent  upon  the  size,  placement  and  material 
charactQvl sties  of  these  structuai  members. 

It  is  recommended  that  adjustments  to  the  inaxiimiiii  surface  temperature  of 
the  reentry  vehicle  be  made  based  on  the  initial  surface  temperatures  provided 
in  Table  8.2.  The  results  presented  in  Section  7 assumes  a 660“R  initial 
surface  temperature,  determine  the  initial  temperature  from  the  intermediate 
results  for  the  particular  region  and  skin  thickness,  Compute  the  difference 
between  this  value  and  660°R.  If  the  initial  temperature  is  less  than  660®R, 
then  subtract  this  difference  from  the  maximum  surface  temperature  determined 
from  Section  7.  Accordingly,  if  the  Initial  temperature  is  greater  than  660‘’R 
then  add  this  difference  to  the  maximum  surface  temperature. 
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APPENDIX  A 

METHODS  AND  PROCEDURES  FOR 
CALCULATING  AERODYNAMIC  COEFFICIENTS 


The  basic  theory  used  In  determining  the  aerodynamic  characteristics  used 
In  this  study  can  be  attributed  to  the  methods  and  procedures  proposed  by 
Allen  (Ref.  A.l).  Jorgensen  (Ref.  A. 2)  used  Allen's  methods  In  computing  the 
normal  force,  axial  force,  and  pitching  moment  coefficients  for  certain  missile 
type  configurations  flying  at  angles  of  attack  between  0 and  180  degrees. 

The  equations  for  determining  the  normal-force,  axial-force,  and  pitching 
moment  coefficients  are  presented  In  Equations  A.  1 - A.  5.  Sketch  (a)  Is 
presented  In  order  \,o  clarify  the  sign  convention  used. 


C » 
m 


and 


Sk(>t,(;l'  (.i) 


L..  - V*  cos*  «' ; 


C*  • C, 


'a  “ Ca, 


180® 


cos*a' ; 


0®  < a < 90° 
90°  < u < 180° 


[V-A.(ll-x,„)1  A„  (Kr^n\ 

J 2 '^^d^  Tp- 


2 . 


0°  < tx  i 90° 


A1 
A. 2 

A. 3 


A. 4 
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sin  2 a'  cos 


+ nC, 


A.,  - x^\ 

i-  [ -T-) 


sin 

90" 


< a 


IBO' 


A. 5 


where 


ex'  ■ a for  0®  < a ;<  90®  and  a'  = 180®  -a  for  90®  5_  180® 

The  aerodynamic  force  center  is  given  by 


ac 


A.6 


The  first  terms  in  liquations  A.i,  A, 4 and  A. 5 come  from  slender-body 
potential  theory.  The  second  terms  represent  the  viscous  crossflow  or  cross- 
flow  attributed  to  flow  separation. 


Crossflow  drag  coefficient  - In  the  expressions  for  and  (equations 

A.I,  A. 4 and  A. 5)  C.  is  the  crossflow  drag  coefficient  for  a section  of  an 

n 

"infinite"  length  cylinder  placed  normal  to  an  uirstream.  It  is  a function  of 
both  the  Mach  number  and  Reynolds  number  components  that  arc  normal  to  the 
cylinder  longitudinal  axis,  and  hence,  for  a body  at  angle  of  attack  it  is  a 
function  of 


and 


sin  o 


A. 7 


Ref^  = Re  sin  a A. 8 

is  commonly  called  the  crossflow  Mach  number  and  Re|^  the  crossflow  Reynolds 
number. 

For  circular  cylinders,  plots  have  been  prepared  for  the  variation  of 
with  M[^  and  Rg|^.  Figure  A.I  gives  the  variation  of  C.  with  M over  the  Mach 
number  range  of  7 or  8 down  to  almost  0. 


A2 


RTI>  0?7-] 


used  to  successfully  predict,  for  most  engineering  purposes,  the  aerodynamic 
characteristics  of  bodies  of  revolution  at  subsonic  Mach  numbers.  For  bodies 
at  supersonic  and  hypersonic  Mach  numbers,  n probably  can  be  assumed  to  bo 
unity,  an  assumption  Indicated  as  being  essentially  correct  from  past  investi- 
gations. 


Axial “force  coefficients  at  angles  of  attack  of  0"  and  180°  - To  pi-edict 
the  variation  of  with  angle  of  attack  by  Equations  and  A. 3,  either 
computed  or  measured  values  of  axial-force  coefficient  at  « « 0“  or  180"  can 
be  expressed  by 


C 


A 


A.  9 


where  Ca  represents  the  wave  or  pressure  contribution  from  the  nose  or  forward- 
facing  base;  C»  is  the  sk1n-fr1ct1on  contribution;  and  C.  Is  the  base- 

Agp  Ajj 

pressure  contribution.  There  are  theoretical  expressions  presented  In  Ref.  A. 2 
that  can  be  used  to  determine  C.  , C.  . C*  . For  this  study  however,  the 

Ay  Agp  Ap 

experimental  data  presented  In  Figs.  A. 3 and  A. 4 were  used  to  determine  the 
axial-force  coefficients  at  0“  to  180°.  These  values  wore  used  In  place  of 
Equations  A. 2 and  A. 3 because  of  the  discrepancy  between  the  theoretical  pre- 
dictions and  the  experimental  measureii^nts. 

For  the  six-degree-of-freedom  trajectory  analysis,  the  aerodynamic  coefficients 
were  computed  for  supersonic  Mach  numbers  only.  Since  the  variation  In  coefficients 
with  Mach  number  in  this  range  was  very  slight,  a constant  coefficient  assump- 
tion allowed  for  many  more  cases  to  be  considered  with  the  trajectory  program. 

This  was  due  to  the  decrease  in  computer  run  time  because  of  not  having  to 
consider  Mach  number  dependent  coefficients.  This  approach  was  entirely  valid 
for  the  heating  analysis  since  any  significant  heating  would  occur  at  tl)e  liii|her 
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Mach  numbers.  The  approach  was  also  valid  for  tho  vehicle  dynamics  at  super- 
sonic Mach  numbers.  When  the  crossflow  Mach  number  (Mn^)  becomes  less  than  one, 
the  approach  starts  becoming  less  valid  since  the  aerodynamic  coefficients  will 
vary  considerably  In  subsonic  flow. 

In  Figs.  A. 3 and  A. 4,  the  computed  values  of  C^,  and  as  a function 
of  angle  of  attack  a are  compared  with  the  experimental  results  for  the  two 
fineness  ratio  cylinder  configurations  and  the  three  fineness  ratio  cone-cylinder 
configurations  shown  below  In  Table  A.l, 


-a:.::.... 


s,- 


Body 

No. 

a 

m 

cF 

V 

F 

^^c 

r 

— 

''s 

d^ 

Nose 

Shape 

1 

6 

3 

6.000 

4.710 

3.000 

19.63 

Flat 

2 

8 

4 

8.000 

6.280 

4.000 

25.92 

Flat 

3 

7 

3.5 

5.500 

3.925 

4.183 

17.31 

Cone 

4 

9 

4.5 

7.500 

5.495 

5.200 

23.62 

Cone 

5 

11 

5.5 

— J 

9.500 

7.065 

6.211 

29.91 

Cone 

Table  A.l 


Bodies  For  Which  The  Aerodynamic  Characteristics  Were  Computed 
Generally,  there  Is  good  agreement  between  the  computed  and  measured  results, 
especially  for  the  variation  of  and  with  a.  The  agreement  between  the 
values  of  were  not  nearly  as  good,  however.  Further  comments  on  this 
disparity  Is  made  In  Ref.  A. 2, 


IT 

I'li 


J 


RTR  027-1 


Body  1 


A/d-8 


Computed 


Measured 

Data 


0 

G 


|l 


Is. 


Fig,  A. 3 Comparison  of  computed  with  measured  aerodynamic  characteris- 
tics for  cylindrical  bodies  with  completely  blunt  noses; 

M - 2.B6. 

m 
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Fig.  A. 4 Comparison  of  computed  with  measured  aerodynamic  character- 
istics for  cylindrical  bodies  with  conical  noses  of  fineness 
ratio  3:  Mo,  ■ 2.86. 
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NOMtNCLATURE 

cross-sectional  (reference)  area  of  cylindrical  portion  of  body 
body  base  area  (at  x -a) 
planforiii  area 

surface  wetted  area  P 
axial-force  coefficient, 

F, 

axial-force  coefficient, 

a P 

crossflov/  drag  coefficient  of  cylinder  section,  jprrn — 

cy^°cy 

pi tching-inoiiient  coefficient  about  station  at  X|j|  from  nose, 
pitching  moment 

normal -force  coefficient,-  p- 

pressure  coefficient, 

body  cross-section  diaiiieter 
body  length 
nose  length 
aftcrsoction  length 

Mach  number  component  normal  to  body  axis,  sin  u 

frce-strodiii  Mach  number 
pressure 

free-streaiii  static  pressure 

dynamic  pressure  component  in  body  axis  diroctinn,  cos"  a 
dynamic  pressure  component  normal  to  body  axis,  g,.,  sin^  k 
free- stream  dynamic  pressure,  ]/?  p 
nose  or  body  cross-section  radius  at  base 

free-stream  Reynolds  number  based  on  body  cross-sect  Ion  diameter, 

QJ-A 

P 

Reynolds  number  component  normal  to  body  axis.  Re  sin  a 
body  volume 
free-stream  velocity 
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X axial  distance  from  body  nose 

X.-  distance  from  nose  to  aerodynamic  force  center 

Civ 

X.  distance  from  nose  to  centroid  of  body  planform  area 

W 

X|^  distance  from  nose  to  pitching-moment  reference  center 

ot  angle  of  attack 

n crossflow  drag  proportionality  factor 

SUBSCRIPTS 

B base 

U wave  or  pressure 
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APPENDIX  B 

6 DOF  TRAJECTORY  PROGRAM  DESCRIPTION 
(Program  REENTR) 

Program  REENTR  is  a versatile  slx-degree-of-freedom  digital  simulation 
providing  a means  of  Investigating  the  reentry  characteristics  of  uncontrolled 
aerodynamic  bodies.  The  program  was  first  developed  for  and  used  In  an  In 
depth  study  of  the  Space  Shuttle  External  Tank  for  NASA.  It  Is  highly  versatile 

in  that  it  employs  several  major  program  options  Including: 

a.  Choice  of  computers  with  minimal  changes 

b.  Machine  plot  option 

c.  Choice  of  sperical,  oblate,  or  Fischer  ellipsoid  earth  models 

d.  Print  format  option 

e.  Random  variable  selection 

In  order  to  streamline  the  program,  Increase  Its  speed  and  still  provide  the 
needed  accuracy,  the  program  was  simplified  by  eliminating  the  option  on  the 
earth  models  and  choosing  a spherical  earth,  and  by  locking  In  the  print  and 
plot  option,  to  reduce  computer  time  In  checking  each  option  for  each  loop. 

In  this  Appendix  a brief  theoretical  overview  of  the  program  will  be 
presented  followed  by  a brief  description  of  the  3 degrees-of-freedom  program 
which  was  derived  from  It. 

13.1  Program  Organization 

Program  REENTR  is.  organized  Into  three  basic  parts:  (1)  Initialization 
of  parameter  (2)  Dynamic  loop  computation,  and  (3)  loop  and  program  derivation. 
Each  part  will  be  presented  in  some  detail  In  the  following  sections. 

B.2  Initialization 


The  program  Is  Initialized  by  first  storing  all  necessary  data.  This 
includes  atmosphere  model,  vehicular  aerodynamic  coefficients,  and  Initial 
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position  state  vector.  The  atmosphere  and  aerodynamic  coefficients  are  normally 
Input  in  tabular  form  and  are  used  by  table  readouts  as  required  throughout 
the  simulation.  Position  data  is  input  for  each  case  and  either  re-1nput  with 
changes  or  re-verified  for  each  subsequent  case. 

The  initialization  portion  is  geared  to  determine  the  state  of  the  vehicle 
at  t ■ 0 in  altitude,  direction,  and  velocity.  This  is  determined  from  the 
following  data  inputs: 

Altitude  (h^) 

Latitude  (Vq) 

Longitude  (X^) 

Relative  Velocity  Magnitude  \ 

Relative  Heading  ® 

Relative  Flight  Path  Angle  (6^^^  ) 

Bank  Angle  ((j)g|^  ) ° 

Aerodynamic  Rol?  Angle  ) 

1 CXq 

Total  Angle  of  Attack  (01,  } 

0 

Since  a spherical  earth  is  assumed,  the  initial  position  vector  is  determined 
as  follows: 


iRo  I “ R + ^ 

where  R is  the  radius  of  the  earth,  then 


^0 


R' 


cos 

cos 

IK  ' 

0 

0 

1 

cos 

T 

sin 

0 

0 

yo 

sin 

4' 

R,  I 

0 

1 i'o 

The  earth  centered  inertial  (I)  frame  Is  thus  established  and  a translational 


vector  to  the  centor-of-grav1ty  of  the  body  detormined. 

A body  centered  reference  frame  is  determined  as  a reference  point  for 
the  velocity  vector.  This  reference  frame  Is  termed  the  geodetic  (D)  frame  and 
1s  defined  as  a reference  frame  with  a north  pointing  horizontal (X)  axis 
with  respect  to  the  Inertial  (I)  frame  and  having  Its  Z-ax1s  pointing  toward 


the  center  of  the  earth.  This  frame  is  always  centered  on  the  body  and  Is 

always  maintained  horizontal  to  tite  surface  of  tlio  earth  and  oriented  as  defined. 
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This  Is  shown  In  Fig.  B.  1 . 

The  direction  cosine  matrix  which  transforms  points  in  the  D frame  to  the 
I frame  Is  determined  from 

A„‘“  = [Xj, 

The  relative  velocity  vector  (R)  frame  is  determined  from  ],  ^^^^1 

Yrgijj*  The  vector  is  determined  by: 

( I'^relol  ’'i-elo  °rel„  | 

'^rslo  “ l'''relo'  ''relo  "relo 

In  order  to  transform  this  to  inertial  velocity,  the  rotation  of  the  earth  must 
be  accounted  for.  Therefore, 


Where  lo^.  is  the  rotational  velocity  of  the  earth,  the  total  Inertial  velocity  vector 
is  then: 


The  orientation  of  the  body  is  determined  by  the  given  body  angles  with  respect 
to  the  velocity  vector.  A body  (B)  to  relative  velocity  vector  (K)  frame  direction 
cosine  matrix  is  determined  from: 

^0^*^  “ 

The  R frame  can  be  related  to  the  D frame  by: 

The  body  to  inertial  transformation  then  becomes; 


A 
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IB 
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With  the  body  orientation  determined,  the  only  other  data  required  Is  the  mass 
characteristics  of  the  reentry  body  Including: 

Mass 

Center-of-gravlty  location  with  respect  to  a known  body  point 
Moment-reference  point  (for  aerodynamics)  from  same  body  point 
Body  dimensions  and  body  rates 

The  Aq  iMtrIx  Is  converted  to  quaternion  (to  be  discussed  In  more  detail  later) 
and  the  dynamic  mon:ent  arm  determined: 


!r^  B 
mrp 


^^mrp"^cg^  ^ '^ref 
^^ffirp'^cg^  ^ '^ref 


The  dynamic  loop  can  now  conmence. 


B.3  Dynamic  Loop  Computations 

As  the  body  reenters  the  atmosphere,  gravitation,  atmosphere,  and  motion 
effects  must  all  be  updated  on  a frequent  basis.  The  dynamic  loop  updates  by 
determining  the  derivatives  of  the  primary  parameters  and  then  Integrating  to 
determine  the  updating  Increment.  The  Integration  method  used  Is  a four  pass 
Runge  Kutler  method.  The  critical  derivatives  are: 

Rotational  acceleration  m 

JitT 

Translational  acceleration  V 

■A  T 

Position  change  R 

Quaternions 

In  the  First  pass  through  the  loop  on  each  ease,  all  of  the  initial  pftrameten. 
are  revorlfled.  In  order  to  simplify  and  compact  the  simulation,  a set  of 
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parameters  known  as  Quaternion,  from  Hamiltonian  calculus,  are  used.  The 
quaternions  are  used  to  represent  the  main  direction  cosine  matrix  and 
are  defined  from  the  Initial  parametersof  that  matrix.  The  calculation  of  these 
parameters  Is  tedious  and  will  not  be  shown  here.  Basical ly  they  are  defined  by: 

Qi  “osIn  U/2) 

Q2  ■ 3s1n  ((|i/2) 

Qj  =YSln  4/2) 

Q4  “ cos  {0  2) 

where  a,  B,  and  y are  elgenaxis  direction  cosines  and  ifi  Is  the  elgenaxis  rotation 
angle.  Integration  and  recomputation  of  components  of  the  matrix  allows 
continual  update  of  vehicle  body  orientation  and  rotation. 

Using  the  above,  the  program  updates  on  a frequent  basis  updating  position, 
velocity,  and  rates  of  each  body  In  contact.  The  forces  on  the  body  are 
determined  based  on  atmospheric  characteristics,  aerodynamic  characteristics,  etc. 
and  force  coefficients  and  moment  coefficient  determined.  The  final  body 
torques  and  forces  are; 

^ ^ref  ^ref 


where 


tnin^  Is  the  torques  on  the  body  based  on  all  forces  supplied. 

A ^ ref  xyz 


where  "£*  Is  the  forces  on  the  body.  The  derivatives  tlieii  are 
xyz 

\T  n r‘  (fA*^-WXlJ) 

I 1s  the  Inertia  matrix  for  the  body 
(j  Is  the  rotational  rate  of  the  body 
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where 


ff  - 1/2 


Q,  -Q,  Q, 

-Q,  Q,  Q, 

/Q»  -Q,-Q, 

I - I 

* -4.  #«  * 


m Is  vehicle  mass 

Is  the  gravitational  vector 

All  derivatives  are  converted  to  Inertial  frame  for  Integration. 


B.4  Loop  and  Program  Termination 

At  the  completion  of  each  Integration  step,  all  parameters  are  updated 
and  the  program  continues  until  stopped.  The  computer  program  checks'  to  de- 
termine print  and  plot  Intervals  and  stores  all  required  data  or  prints  it. 
At  program  termination,  the  final  data  Is  printed  and  plotted  as  required. 


B.5  Three  Degree-Of-Freodom  Program 

A 3 degree-of-freodoiii  program  was  also  developed  for  the  Shuttle  studies 
and  used  1n  this  study.  This  program  Is  an  adaptation  of  the  six  deqree-of- 
freedom  program  employing  a point  mass  lift/drag  model  based  on  vehicle  aero- 
dynamics. 

The  primary  advantage  to  the  use  of  this  model  Is  its  high  speed  (Approx. 

1 second  per  100  seconds  of  trajectory)  thus  allowing  the  I’un  of  large  numbers 
of  similar  cases,  perhaps  varying  one  parameter  over  a large  range.  Where  the 
reentry  bodies  trim,  this  becomes  a high  speed  way  of  studying  the  changes 
in  reentry  characteristics  with  eight  variations  In  atmospheric  conditions,  etc. 


BIO 


B.5  CONCLUSION 


The  two  programs  described  In  this  Appendix  wore  modified  slightly  to 
provide  data  useful  In  computing  aerodynamic  heating  data.  Direct  type  output 
of  many  trajectories  could  be  fed  directly  Into  heating  calculation  programs 
to  determine  heating  characteristics.  Program  RECNTR  and  Its  derivatives  are 
currently  available  In  REMTECH*s  files  for  future  studies  of  this  type. 
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MATHEMATICAL  NOTATION 


NOTATION 

vl 


DEFINITION 

Vector  resolved  into  I-frame  coordinate  axes.  Axes  systems 
used  arc: 

1 - Earth  centered  inertial 
D - Local  normal,  north  pointing 
G " Geocentric  local  normal,  north  pointing 
R - Relative  velocity  vector,  local  horizontal 
B - Body  fixed 

Direction  cosine  matrix  which  transforms  vectors  from  the  y-fratne 
to  the  X- frame 

Initial  value  of  X 


dx/dt 

1 0 0 " 

0 cosO  -sino 
0 sInG  cos  0 

: H 

COSO  0 sinO 

0 10' 
■sInO  0 cosG 

COSO  -s1nO  0 

sInO  COSO  0 

0 0 1 


Direction  cosine  matrix  for  x-oxis  rotation 
through  angle  0 

Direction  cosine  matrix  for  y-ax1s  rotation 
through  angle  0 

Direction  cosine  matrix  for  z-axis  rotation 
through  angle  0 


NOTE;  All  direction  cosine  matrices  used  herein  transform  from  the  new 
coordinate  frame  to  the  original  coordinate  frame. 
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!:■ 

LIST  OF  SYMBOLS  T. 

■Ii: 


SYMBOL 

COMPUTLR  SYMBOL 

PFFINITION 

S 

CA 

Aerodynniiilc  axial  fores  coefficient 

s„ 

CM 

Aerodynamic  pitching  mment  coefficient 

s 

CN 

Aerodynamic  normal  force  coefficient 

^llp 

CLP 

Rolling  tiMimont  damping  coofficlent 

^lliq 

CMQ 

““  Pitching  moment  damping  coefficient 

4 

^nr 

CNR 

Yawing  moment  damping  coefficient 

I'niriy 

CLMNT 

Total  aerodynamic  moment  coefficient 
vector  with  respect  to  c.g. 

•(< 

CPLMN 

Total  aerodynamic  moment  coefficient 
vector  less  c.g.  olfset  offect 

CXY2 

Total  aerodynamic  force  coefficient  vector 

^rof 

ORL'- 

Vehicle  diameter 

1*  15 
'a 

FAB 

Aerodynamic  force  vector 

i'h 

g 

C.MAO 

Magnitude  of  acceleration-of-gravi ty 

g 

G 

Accoleratlon-of ‘gravity  vector 

OM 

GM 

Larth  gravitational  constant 

h 

H 

Altitude 

I 

IXY/. 

Vehicle  moment  of  inertia 

• 1 

1 

IXY/I 

Inverse  moment  of  Inertia 

^■of 

LRir 

Vehicle  length 

III 

GMASS 

Vehicle  mass 

M.\ch 

MACH 

Mach  lumiher  (V,.j,^/Vs) 

P 

pqR(l) 

Roll  rate 

q 

pqR(;>) 

Pitch  rale  (■=^11^) 
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SYMBOL 

COMPUTER  SYMBOL 

r 

PQR(3) 

q 

QBAR 

Q iQ  iQ  >Q 

2 » •* 

Q(1).Q(2).Q(3),Q(4) 

R 

R 

>*5 

RS 

Jjef 

SREF 

T® 

A . 

TAB 

Iv  , 1 

/rel 

VRMAG 

* rel 

VRELX 

h 

vs 

^x 

w 

vwx 

^cg’^cg’^cg 

XCG.YCG,ZCG 

X Y Z 

''mrp,  mrp,  mrp 

XMRP.YMRP.ZMRP 

“T 

ALPHT 

"'^rel 

GAMR 

\ 

LAMDA 

P 

RHO 

°rel 

SIGR 

01 

PHIA 

‘’’bk 

PHIBK 

^D 

PSID 

OMGE 

0) 

PQR(I) 

DEFINITION 
Yaw  rate  (-w^) 

Dynamic  pressure 
Quaternion  parameters 
Earth's  radius 

Kennedy  Space  Center  earth  radius 
Vehicle  reference  area 
Aerodynamic  torque  vector 
Magnitude  of  relative  velocity  vector 
Relative  velocity  vector  In  X- frame 
Speed  of  sound 

Wind  velocity  vector  In  x-frame 
Reference  location  of  c.g. 

Moment  reference  point  location  with 
respect  to  c.g. 

Total  angle-of-attack 

Relative  flight  path  angle 

Longitude 

Atmospheric  density 
Relative  heading 
Aerodynamic  roll  angle 

Bank  angle 
Geodetic  latitude 
Earth  rotational  velocity 
Vehicle  rotational  rates 
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APPENDIX  C 

REENTRY  HEATING  METHODOLOGY 

The  aerodynamic  heating  methods  used  In  this  study  were  those  contained  in 
the  GENTRY  computer  code  (Ref.  C.l)  and  modifications  for  this  specific  class 
of  problems.  A discussion  and  definition  of  the  aerodynamic  heating  methods 
used  for  each  part  of  the  payload  geometry  is  given  in  this  appendix.  The 
geometries  considered  were  flat  face  right  circular  cylinders  and  cone-cylindor- 
flat  face  bodies.  The  heating  methods  for  each  geometry  is  described. 

Free  stream  properties  were  obtained  for  each  trajectory  point  using  the  1962 
standard  atmosphere.  All  thermodynamic  and  transport  properties  of  air  were 
obtained  using  the  method  of  Hansen  (Ref.  C.2).  Equilibrium  thermodynamic 
properties  of  air  which  include  its  major  components  (Oa*  Na.  0,  N,  0'^,  N'*',  e“) 
were  obtained  in  closed  form  from  a set  of  approximate  partition  functions. 

The  partition  functions  and  thermodynamic  relations  were  programed  such  that 
the  total  enthalpy,  entropy,  speed  of  sound,  average  molecular  weight,  heat 
capacity  and  species  concentration  can  be  found  for  a given  temperature  and 
pressure.  The  calculation  of  transport  properties  by  Hansen's  method  -is  based 
on  simple  kinetic  theory  of  hard  spheres.  The  viscosity  is  a function  of 
species  density,  mean  velocity  and  mean  free  path.  The  thermal  conductivity 
is  computed  accounting  for  energy  transfer  by  molecular  collision  and  for  energy 
transfer  by  specie  diffusion. 

A boundary  layer  transition  criterion  which  is  simple  and  appropriate  for 
design  work  was  used.  The  criterion  can  be  stated  as: 
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I Laminar  If  Ret„j|  < 100000 

Transition  Transitional  If  100000  < < 200000 

Criterion 

Turbulent  If  Re^^j  > 200000 

The  atmospheric  properties,  air  thermodynamic  and  transport  properties  and 
transition  criterion  were  applied  to  all  geometries. 

CYLINDER  GEOMETRY 

The  specific  calculations  procedure  for  the  cylinder  Is  determined  by  the 
angle  of  attack  and  the  transition  criterion.  The  following  bar  chart  gives 
reference  to  the  heating  methods  used. 


Laminar 

a “ 

Turbulent 


Fig.  C.l  Stagnation  Line  Cylinder  Heating  Methods 
The  numbers  In  Fig.  C.l  refer  to  tables  which  define  the  heating  method. 
All  the  methods  require  a recovery  enthalpy  calculation.  The  recovery 
factor  Is  computed  using  the  Newtonian  type  theory  presented  In  Table  C.7. 

In  order  to  utilize  these  heating  methods  the  local  flow  conditions  are 
required.  The  following  steps  are  used  to  compute  tlie  flowfleld. 

Zero  Angle  of  Attack  Flow  Fields 

1.  Normal  shock  calculation  Is  made'. 

2.  Total  post  shock  conditions  ore  computed. 

3.  Using  the  stagnation  entropy,  the  flow  Is  expanded  to  free  stream 
pressure. 

4.  Edge  conditions  are  determined  from  the  stagnation  entropy  and 
edge  pressure. 

C2 
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Angle  of  Attack  Flow  Fields 

If  the  angle  of  attack  plus  the  body  angle  Is  greater  than  1.0  degrees 
an  angle  of  attack  calculation  is  made  using  the  following  steps: 

1.  An  oblique  shock  calculation  Is  made  for  a shock  parallel  to  the 
local  body  angle  If  the  normal  component  of  free  stream  Mach 
number  Is  greater  than  one. 

2.  If  the  normal  component  of  Mach  number  Is  less  than  one.  the  edge 
enthalpy  and  pressure  are  computed  by 

Edge  enthalpy  hg  « h^  + U»*^/2gJ 

Edge  Pressure  Pg  “ Pa,  (1  + 

where  the  subscript  n means  normal  component.  Based  on  these  two 
properties  the  remaining  edge  conditions  are  computed. 

3.  If  an  oblique  shock  calculation  was  made  (step  1).  the  edge  conditions 
are  computed  by  Isentropically  compressing  the  normal  component  of 
velocity 

Edge  enthalpy  Sg  » 

Edge  enthalpy  hg  ■ h^  + U^*j,/2g0 

6 ■ post  shock 

Typical  results  using  the  preceding  heating  methods  and  flow  field  methods  are 
compared  with  experimental  data  In  Fig.  C.l. 

Heating  Distributions 

In  order  to  compute  the  heat  transfer  coefficient  distribution  for  a 

tumbling  or  rolling  case  a method  based  on  data  was  used.  This  distribution  1s 

given  In  Table  C.8.  Thus  with  the  stagnation  line  value  from  the  previous 
methods  and  the  local  to  stagnation  line  heat  transfer  coefficient  ratio,  the 
local  value  can  be  obtained.  The  heating  is  then  computed  using 
q “ h(H^-H^^) 

where 

q = local  heating  rate 
h » local  heat  transfer  coefficient 
“ local  recovery  enthalpy 
11^^  » wall  enthalpy 
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The  base  geometry  was  approximated  as  a flat  face  of  a right  circular 
cylinder.  The  specific  calculation  procedure  for  the  base  Is  determined  by  the 
angle  of  attack  and  transition  criterion  and  Is  shown  pictorically  below: 


Turbulent 


Fig.  C.3  Base  Flow  Heating  Methods 

The  flow  field  calculations  for  the  different  heating  methods  were  made  as 


described  below: 


a < 95**  and  g > 175° 


1.  Normal  shock  calculation  Is  made. 

2.  Total  post  shock  calculations  are  computed  and  used  as  edge  conditions 


Flow  Fields  for  9S°  < a < 175^ 

An  oblique  shock  calculation  Is  made  for  a shock  parallel  to  the  base 
surface  If  the  normal  component  of  free  steam  Mach  number  Is  greater 
than  1.0. 

If  the  normal  component  of  Mach  number  Is  loss  than  one,  the  edge 
enthalpy  and  pressure  are  computed  by 

Edge  enthalpy  hg  = + U„%/2cjJ 

Edge  pressure  Pg  * Pog  (1  + .211/, 

where  the  subscript  n means  normal  component.  Based  on  these  two 
properties  the  remaining  edge  conditions  are  computed. 

If  an  oblique  sho''k  calculation  was  made  (step  1),  the  edge  conditions 
are  computed  by  Isentropically  compressing  the  normal  component  of 
velocity 
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The  heating  1$  computed  to  the  center  of  the  base  and  assumed  to  be  constant 
over  the  entire  base. 

CONE  GEOMETRY 

A separate  cone  calculation  methodology  was  not  required  since  the  cone 

heating  may  be  related  to  cylinder  calculations.  The  Stanton  number  data  and 

theory  presented  by  Ref.C.3i  shown  In  Fig.  3.  demonstrates  that  swept  cylinder 

heating  theory  can  be  used  for  cone  heating  calculations.  The  heating  rate 

for  a cone  related  to  the  cylinder  Is  simply  a function  of  the  radius  ratio. 

qcone  -/RcyL.\« 
qcyl  \ Rcone  / 

a " 0.5  laminar 

a ■ 0.2  Turbulent 

The  local  cone  radius  Is  used. 

All  of  the  heating  calculations  for  this  report  were  made  using  reference 
conditions.  These  are 

Surface  Temperature  » 660* R 
Cylinder  Diameter  ■ 24  inches 

Precedures  are  described  In  the  main  text  to  adjust  for  diameter  and  surface 
temperature  effects  on  heating  rates  and  load. 
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Table  C.1 

ECKERT'S  LAMINAR  FLAT  PLATE  HEATI4J6  METHOD 
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Eckert's  laminar  flat  plate  heating  relation  can  be  expressed  as: 
h » 0.332  9c 

\ri^  (lbm/ft*sec) 

where 

Dw  a e--MeX  Reynolds  no.  evaluated  at 

e v*  reference  conditions 

The  reference  conditions  are  obtained  from  the  reference  enthalpy  [i 

H*  - He  + 0.5  (Hw  - Hg)  + 0.22(Haw  - Hg) 
and  edge  pressure.  Other  reference  properties  are  evaluated  as: 

P*  - f (H^  Pe) 

11*  “ f (H*.  Pfl) 


For  geometries  other  than  a flate  plate,  the  Mangier  transformation 
may  be  evaluated  as: 

Tm  “ 1 Flate  Plata 

= 2 Cylinder 

" 3 Cone 

Nomenclature 

Qj  *=  32.174  Ibm/slug 

h » Heat  Transfer  Coefficient  (Ibm/ftHcc) 

Hp  « Edge  Enthalpy  (Btu/lbm) 

Haj^  » Adiabatic  Wall  Enthalpy  (Btu/lbm) 

“ Wall  Enthalpy  (Btu/lbm) 

Pg  = Edge  Pressure  (atm.) 

Ug  = Edge  Velocity  (ft/sec) 

X " Surface  Distance  From  Origin  to  Point  of  Interest 
p*  “ Reference  Density 
M*  a Reference  Viscosity 
‘m  “ Mangier  Transformation 
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Table  C.2 

STAGNATION  POINT  HEATING 
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Stagnation  point  heating  Is  based  on  the  theory  of  Fay  and  Riddell 
(Ref.Tl). 


where 


9 “ ^ 


W (Ht“  V 

W fu 


0.76 


0.4 


1 + 


(Le'*-  1)  Hi 


a sphere  and  1,22 
a cylinder  and  M^>  1 .12 


» U (1.5  - 0.378M'‘  - 0.02625  M**  ) for  a sphere  and  M < 1.22 

_JB  00  60  I" 

R 

- (2.0  - 0.87, ’M'  - 0.328M»  ) for  « cvllndsf  and  t|,  < l.)2 

and  the  dissociation  enthalpy 

“d  ■ 

h°  • 6636.26  Btu/Lbm 

h“  - 14456.53  Btu/Lbm 
N 


Nomenclature 


q ■ Heating  Rate  (Btu/ft'^sec. ) 

“ 1 (slug  ft./lb^sec.*) 

g'  « 32.174  (Lbm/slug) 

" Density  (slug/ft.’) 

U " Viscosity  (LbrSec/ft'^) 

U ■ Velocity  (ft. /sec.) 

P ■ Pressure  (Lb-/ft.O 

H ■ “ Enthalpy  (Btli/lbni) 
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Table  C. 2 (Cont.) 
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Lg  = Lewis  Number 

= Dissociation  enthalpy  of  air  (Rtu/Lbin) 

" Mass  fraction  of  species  1 (1  ■ 0,  oxygen  and  1 » N,  nitrogen) 

R “ Nose  radius  (ft.) 

■ Prandtl  Number 

Subscripts 

e ■ Edge 

t ■ Stagnation  (post  shock) 

w ■ Wall 

00  ■ Fi^ee  stream 


References 

T1 . Fay,  J.  A.  and  Riddell,  F.  R.,  "Theory  of  Stagnation  Point  Heat  Transfer 
In  Dissociated  Air",  J.  Aeronaut.  Sci.,  Vol . 25,  No.  2.  Feb.  195P.  dp. 
78-85,  121. 
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Table  C.3 

LAMINAR  STAGNATION  LINE  HEATING  TO  YAWED  CYLINDERS 


The  laminar  stagnation  line  heating  by  the  method  of  Sayano  and  Greenwald 
(Ref.  T1)  Is  expressed  by  the  following  empirical  relation. 

cos*'*  A + (qfp/q^-o)  s<n*A 

where 

^A-o  “ ^FR 

Nomenclature 

h ■ Heat  Transfer  Coefficient 

hpp  « Fay  and  Riddell  Spherical  Stagnation  Point 
Heat  Transfer  Coefficient 

q^p  ■ Flate  Plate  Heating  Rate 

q ■ Cylinder  Heating  Rate  at  Zero  Sweep 

A “ Sweep  Angle 

Reference 

T1  Sayano,  S.  and  6.  F.  Greenwald,  "Approximate  Method  for  Calculating  Heat 
Transfer  to  Yawed  Cylinders  1n  Laminar  Flow",  JSR  Vol.  10,  No.  2,  Feb.  1973. 


C9 


RTR  027-1 


Table  C.4 

SPALDING-CHI  METHOD  MODIFIED  FOR  REAL  GAS  HEAT  TRANSFER 
FOR  TURBULENT  BOUNDARY  LAYER  FLOW 


(1)  Spalding-Chl  define  Fg,  F^^.  Fr^  (Ref.Tl)  which  are  functions  of 
Mach  Number  and  Temperature  alone  such  that 

1/2  C/,  • 'l'o(Vte’  ' ♦s"'i<s"Es> 


Where 


Z.(T)  -f  (H.P^)  & H • + 


(3) 


^Rof.Ta) 


(5) 


Where  g(l)  » 9.2808635 
q(2)  ■ -4.7340248 
g(3)  » 6.6858663  10"* 
g(4)  “ -4.1876614  10"* 
g(6)  - -5.5054577  10’" 
g(6)  » 2.83C/291  10"" 

g(7)  - -2.1249608  10"* 
g(8)  - 8.0162000  10"’ 
g(9)  « -1.5900985  10‘* 
g(10)  ■ 1.3236350  10"*® 


Rp  ■ ie^ 
^s  Hg 


» Mangier  Transformation  (Ref.T3) 
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Nomenclature 


fiTURB 

H 

P 

Pf 

ReO.  Re^ 


s 


St 


T 


u 

Z 


P 


Subscripts 

aw 

e 

1 

w 


Skin  Friction  Coofficient 

Heat  Transfer  Coefficient  For  Turbulent  Flow 

Enthalpy 
Pressure 
Prandll  Number 

Local  Reynolds  Number  Based  on  Momentum 
Thickness  0 and  Characteristic  Length  s 
Respectively 

Inverse  of  Von  Karman  Reynolds  Analogy  Factor 

Stanton  Number 

Temperature 

Velocity 

Compressibility 

Density 


Recovery 

Local 

Incompressible 

Wall 


References 


Tl.  Spaldingi  O.B.i  and  Chi,  S.W.,  “The  Drag  of  A Compressible  Turbulent  Boundary 
Layer  on  a Smooth  Flat  Plate  With  and  Without  Heat  Transfer",  Journal  of 
Fluid  Mechanics,  Vol.  18,  Pert  I,  pp,  117-143,  Jan.  1964. 

T2.  Wallace,  J.E.,  "Hypersonic  Turbulent  Boundary  Layer  Studies  at  Cold  Wall 
Conditions",  1967  Heat  Transfer  and  Fluid  Mechanics  Institute,  La  Jolla,  CA 
June  1967. 

T3.  Komar,  J.J.,  "Improved  Turbulent  Skin-Friction  Coefficient  Predictions  Util- 
izing the  Spalding-Chl  Method",  Douglas  Aircraft  Company,  Douglas  Report 
DAC-59801,  Nov.  1966. 
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Table  C.5 

TURBULENT  STAGNATION  LINE  HEATING 
NEAR  a « 0 


Beckwith  and  Gallagher's  stagnation  line  theory  predicts  zero  heat  transfer 
coefficients  at  an  angle  of  attack  of  zero.  In  order  to  achieve  a transition  to 
more  realistic  value  the  following  approximation  for  turbulent  flow  Is  used. 

0 ^ o 1 5 deg. 

0 ■ 18  (5  - o) 

“ hgQ  4 h^  . 0 sin  0 

Nomenclature 

h - Heat  transfer  coefficient 

h „ - Heat  transfer  coefficient  at  a ■ 0,  computed  using  flat 
plate  conditions  processed  by  a normal  shock  and  expanded 
to  edge  pressures. 

a - Angle  of  attack 

hgg  - Beckwith  and  Gallagher  heat  transfer  coefficient 
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BECKWITH  AND  GALLAGHER  TURBULENT  REAL  GAS 
YAWED  CYLINDER  STAGNATION  LINE  HEATING 


g 

h 

p 

Pr 

R 

u 

X 

A 

p 


32.174 

Lbm/slug 

Heat  Transfer  Coefficient 

Lbm/ft*sec 

Pressure 

Lbf/ft® 

Prandtl  Number  edge  conditions 

Dimensionless 

Radius  of  Cylinder 

Feet 

Velocity 

Ft/sec 

Distance  from  Stagnation  Line 

Feet 

Sweep  Angle 

Degrees 

Viscosity 

Lbm/ft-sec 

Density 

Clugs/ft’ 

Subscripts: 


SL  St-icjnation  Line 
0 Edge  Conditions 
^ Total 


Superscript- 

* Evaluate  at  Eckert's  Reference  Enthalpy  & edno  Pressure 


Reference: 

1.  Beckwith,  1.  E.,  and  Gallagher,  J.  0.,  "Local  Heat  Transfer  and  Recovery 
Teniperaturos  on  a Yawed  Cylinder  at  a Mach  Nuiiibei  of  4. IB  and  Hitjh  Rey- 
nolds Number,  NASA  TR  R-104,  1961. 
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TABLE  C.7 

RECOVERY  FACTOR  DISTRIBUTION 
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The  local  recovery  factor  or  recovery  enthalpy  Is  computed  using  the 
following  Newtonian  type  of  distribution  (Ref.  Tl). 


hj,  ■ h„  + [r^  sin*  a cos*  0 + r^^,Q  (sin*  a sin*  0 f cos*  a)]U4/2gJ 


where 


r ■ 


*^0  ■ 0 


U*/2gJ^ 


U^2gJ 


a ■ 0 


1 .0  for  0 ^ 0 <,  90 
,2 


(H^  - A^/2gJ  - h^)  + P*/n  A^/2gJ 
U^2g0 


A|^  ■ Leeside  speed  of  sound  based  on  total  temperature 
n " 2 for  laminar  flow 
n « 3 for  turbulent  flow 
For  subsonic  flow  r^  " 1.0  for  all  0. 


Nomenclature 

hoc  ■ Free  stream  static  enthalpy  (Btu/lbm) 
h>«  " Recovery  enthalpy  (Btu/lbm) 

Hx  ■ Total  enthalpy  (Btu/lbm) 

U„  • Free  stream  velocity  (ft. /sec.) 
g ■ 32.174  Ibm  ft./lbf  sec 
J - 778.66  ft.  Ibf/Btu 
r ■ Recovery  factor  0.85  laminar  flow 

0.897  turbulent  flow 
a Angle  of  attack  (degrees) 

0 “ Angular  position  from  stagnation  line (degrees) 
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TABLE  C.7  (cont.) 


The  following  figures  show  a comparison  of  data  and  the  preceding  equation. 


• SH-12F  SRB 
Recovery  Factor 
Data,  clean  skin 
cylinder  values. 


References 

11  C.  D.  Engel,  "Tumbling  Entry  Heating  .‘'Yograms  - Teheat  and  Bentry  Users' 
Manual,"  RTR  OllC-1,  NAS8-21810.  July  1976. 
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Table  c.8 

TURBULENT  CIRCUMFERENTIAL  HEATING 
DISTRIBUTION  AROUND  A CYLINDER 


The  circumferential  heating  distribution  around  a right  circular  cylinder 
has  been  correlated  (Ref.  T1)  using  the  following  equations 

V^sl  ’ ^*0  cos  A + cos  2A 

+ (2Ki  + Ka  - 3)/4 

where 

Kj  ■ © "‘a  cosh*  (3.98  «) 

Ka  ■ © •*“  cosh  (a)  for  0 < a < 0.785398 

k,  ■ © •■»«  cosh  (a)  tl  + 11. 8/(1  " 1.309)j] 

for  0.785398  < a < ii/2 

hy  ■ Undisturbed  heat  transfer  coefficient 

hji  ■ Stagnation  line  heat  transfer  coefficient 

0 ■ Angle  of  attack  resolved  Into  90  degrees  (radians) 

A ■ Circumferential  angle,  from  stagnation  line  resolved 
Into  180  degrees 

A comparison  of  the  correlation  with  data  Is  given  In  the  following  figure. 
References 

Tl.  C.  0.  Engel,  J.  V,  McAnally,  B.  M.  Levine,  "SRB  Data  Book  Procedure 
Documentation",  RM  020-2,  NAS8-31360,  March  1976. 


CIS 


r 


RTR  027- ] 


Table  C.a  (Cont.) 


° 20  i<0  60  80  100  120  n*0  'l60  180~ 

^ (Deg . ) 


i‘i  i 


•i 

i 


Comparison  of  Correlation  Equation  and  SH-I2F 
Test  Data  for  Undisturbed  Heating  (M  ■ 3,75. 
Re»/ft  - 3.5  X I0\  D - I*. 38  In).  “ 
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TABLE  C.9 
LEEWARD  BASE  FLOW 


The  heating  to  the  base  of  a cylinder  when  the  base  fare  Is  leeward  1s 
approximated  by: 

q ■ h (Haw-Hw) 

where 

h ■ 0>05  hpi^ 

hpi^  ■ Hemispherical  stagnation  point  heating 

coefficient  for  cylinder  radius  (Table  C.2) 

Hftw  ■ 0.86  Ht 

Hm  • Wall  enthalpy 
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TABLE  C.10 
STAGNATION  BASE  FLOW 


The  heating  to  the  flat  face  of  right  circular  cylinder  Is  determined 
from  the  theory  of  Fay  and  Riddell  together  with  a velocity  gradient  modification 
from  Bolson  and  Curtiss  (Ref.  T.l).  The  heat  transfer  coefficient  Is  expressed  as 

h - hpR  (C.46628)  M» 

where 

hcD  ■ Hemispherical  stagnation  point  heating  coefficient 
for  the  cylinder  radius  (Table  C.2) 

Noe  ■ Free  stream  Mach  number 

Reference 

Tl.  J.  C.  Bolsoni  H.  A.  Curtiss,  "An  Experimental  Investigation  of  Blunt  Body 
Stagnation  Point  Velocity  Gradient,"  ARS  Journal,  pp,  130-135,  February  1959. 


Heat  Transfer  Coef./Ref.  Heat  Transfer  Coef.(hs/hri0f) 
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APPENDIX  D 

CONDUCTION  HEAT  TRANSFER  ANALYSIS 


The  first  approach  to  predicting  the  surface  temperature  for  the  cases 
of  Interest  was  to  use  a heat  conduction  program  which  accounted  for  all 
possible  circumstances  such  as  variable  material  propert1es»  radiation  at 
the  Interior  and  exterior  surfacesi  transient  effects,  and  variable  convective 
heat  transfer  coefficients.  Reference  D.l  documents  such  a program  which  uses 
a finite  difference  scheme  to  solve  the  heat  transfer  equations: 


ha  (tRa-te)  " (te**-taa‘»)  - ■ (exterior  surface)  (01) 

> ^lc|^y  (steady- state)  (Internal  point)  (D2) 

pc  (transient)  (Internal  point)  (03) 

hi  (tRi-t^)  - e^t  (ti^-t^^**)  - - (k^y  (Interior  surface)  (D4) 


where 

h - local  convective  heat  transfer  coefficient 
to  - recovery  temperature 
t • local  material  temperature 
e - emissivlty 

0 - the  Stefan -Boltzmann  constant 

ta  - temperature  to  which  the  point  radiates 

k - conductivity' 

X - Internal  distance  from  the  surface 

p - density 

c - specific  heat 

e - time 

'jiQ  - Interior  and  exterior  surface 

The  skin  Is  approximated  by  an  Infinite  plate  and  equations  D1  - D4  have  only 
one  space  dimension  -x. 

Preliminary  efforts  dealt  with  choosing  finite  Increments  of  space  and  time 
which  would  give  adequate  solutions  without  using  computer  time  unnecessarily. 

It  was  found  that  use  of  four  Internal  nodes  along  with  the  Interior  and  exterior 
node  gave  results  almost  Identical  to  those  achieved  using  6 or  more  Internal 
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nodes.  Also  time  steps  of  one  second  produced  siisIlAr  results  when  compared 
to  more  rigorous  time  steps  of  0.5  and  0.1  seconds. 

After  several  test  cases  were  run,  It  was  found  that  the  emissivlty  of 
aluminum,  oQ.lSi  was  low  enough  that  radiation  losses  were  very  small.  It  was 
also  found  that  the  temperature  gradient  through  the  skin  was  very  small. 

These  facts  led  to  a considerably  simplified  equation  which  assumed  no  radiative 
heat  transfer  and  no  Internal  resistance  to  heat  flow  (high  conductivity): 

q " cpV^  - hgAe  (tR-t)  (06) 

where 

q - heat  transfer 
V - volume 
A - area 

For  the  one-dimensional  case,  V/A  Is  equal  to  Ad/A  or  d,  the  thickness. 

Solution  of  the  above  equation  yields: 

Nu  ■ e ■ to  ■ temperature  9 e • 0 (06) 

or 

t ■ e"^^®  (ti-tp)  + tR  (D7) 

This  equation  Is  readily  solved  using  a calculator. 

To  compare  the  two  methods  of  solution  three  cases  (Table  Dl)  were  run 
using  1/4  Inch  thick  aluminum  skin.  Each  method  assumes  that  h and  tq  are 
constant  during  the  time  Increment  Involved.  The  second  method  calculated  the 
average  h andtq  for  the  time  step  and  uses  these  values.  The  first  method 
integrates  h andtq  across  the  step.  A listing  of  the  program  used  to  solve 
equation  07  with  a HP-97  calculator  Is  shown  In  Table  D2. 


02 


RTR  027-1 


The  results  are  compared  in  Fig.  Dl.  The  plots  for  all  practical  purposes 
coincide  proving  that  the  assumption  of  negligible  Internal  resistance  was 
correct  for  1/4  Inch  (and  therefore  thinner)  skins. 

The  effects  of  thickness  on  the  accuracy  of  the  second  method  can  be 
seen  In  Fig.  02  where  both  methods  were  used  to  predict  the  surface  temperature 
for  the  worst  case  (Case  3)  when  the  skin  thickness  Is  3/8  inch.  The  second 
method  does  not  respond  as  quickly  as  the  more  correct  first  method,  and  It  Is 
anticipated  that  thicker  skins  would  have  to  be  analyzed  using  the  time- 
consuming  first  method. 

The  effects  of  thickness  on  the  peak  surface  temperature  for  the  three 
cases  are  shown  In  Fig.  03.  As  would  be  expected,  the  thicker  skins  provide 
additional  heat  sinking  so  the  peak  temperature  Is  significantly  lowered. 

In  the  second  method  of  calculation  It  was  assumed  that  the  response  of 
the  skin  to  each  Input  of  heat  during  the  time  step  would  be  slow  enough  that 
the  skin  would  not  come  to  steady  state  before  a new  time  Increment  and  thus 
a new  heat  transfer  rate  began.  To  check  this,  an  extremely  simplified  third 
method  was  Investigated  which  assumed  steady  state  response  at  each  time  step 
using  the  equation: 

q ■ cpV  (tj-to)  (08) 

where,  tg  Is  the  skin  temperature  resulting  from  the  cumulative  heat  transferred, 
q.  to  the  skin  volume.  V.  The  skin  volume  would  be  Aa  where  A Is  equal  to  unity 
for  the  one-dimensional  assumption  made  throughout  the  analysis.  Therefore  the 
temperature  at  anytime  Is: 

‘s  ■ *0  ♦ (D9) 
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To  compare  the  first  and  third  method,  values  of  q were  obtained  using 
BENTRV  (Ref.  C.1, Appendix  C)  which  assumed  a constant  660S'  surface  temperature. 
These  values  and  values  of  c and  p at  200'*C  (IIIB^R)  were  used  to  calculate 
the  peak  temperature  for  the  sample  cases.  The  results  are  shown  In  Fig.  04 
and  Indicate  an  almost  exact  correspondence.  Considering  that  the  BENTRY  q 
calculation  assumed  a constant  temperature  It  was  concluded  that  q had  little 
dependence  on  surface  temperature.  It  was  also  concluded  that  transient 
conduction  effects  were  negligible. 

The  surface  temperature  calculation  scheme  finally  used  was  a form  of  the 
third  method.  First,  using  the  temperature  at  the  beginning  of  the  trajectory 
time  step,  q was  calculated  for  each  step.  Second,  the  temperature  at  the  end 
of  each  time  step  (which  was  also  the  beginning  of  the  next  step)  was  calculated 
using  equation  09.  This  process  continued  until  the  skin  began  to  cool  (negative 
q)  - the  last  calculated  skin  temperature  was  therefore  the  peak.  After  several 
trials  were  made,  It  was  found  that  a variable  specific  heat  was  needed  for 
accuracy  and  the  equation: 

c - 0.08697  + 3.686  X lO"^  - 2.803  X 10"^t*  + 8.23  X 10'^‘t*  (010) 

evaluated  at  the  beginning  of  each  time  step,  was  used  in  all  subsequent 
calculations. 

The  simple  method  of  calculation  used  Is  Justified  by  comparison  with  the 
rigorous  method.  The  method  compares  well  because  thin  aluminum  skins  react 
quickly  to  the  Input  heat  transfer  rates.  Use  of  materials  with  significantly 
lower  conductivities  or  thicker  skins  would  Invalidate  the  assumptions  Implicit 
In  this  method,  and  a more  sophisticated  approach  would  have  to  be  used. 
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TABLE  D2 

PROGRAM  FOR  CALCULATING  SOLUTIONS  ON  HP-97  CALCULATOR 
(Negligible  Internal  Resistance  Assumed) 

Initialization 

StoTe  A0' : — 

Store  hi  — 

Store  tj^  

Computation 


1 

cpT 


Enter  with  ti  (J«1}tj-t.)*  » 

Recall  

Enter  Enter  t^j^^ * 

Average  h for  time  step,  


NuO 

g-Nu® 

Recall  tpj 


Average  tRfor 


time  step. 


*Rj+tRj+l 
— 2 


■Nu0‘ 


Vl  ■ (tj-tple  + t 
Go  to  label  B with  t 


J “J+l 


Finish 


* Temp  at  beginning  of  time  step 
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FIRST  NETHOO 
SEOHiD  METHOD 


nethod  calculations  for  1/4  Inch  aluminum  skin 


Confjarison  of  first  and  second  method  calculations  for  3/8  inch  aluminun 


